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The Road From Classical to Quantum Codes: 

A Hashing Bound Approaching Design Procedure 

Zunaira Babar, Panagiotis Botsinis, Dimitrios Alanis, Soon Xin Ng and Lajos Hanzo 


Abstract —Powerful Quantum Error Correction Codes 
(QECCs) are required for stabilizing and protecting fragile 
qubits against the undesirable effects of quantum decoherence. 
Similar to classical codes, hashing bound approaching QECCs 
may be designed by exploiting a concatenated code structure, 
which invokes iterative decoding. Therefore, in this paper we 
provide an extensive step-by-step tutorial for designing Extrinsic 
Information Transfer (EXIT) chart aided concatenated quantum 
codes based on the underlying quantum-to-classical isomorphism. 
These design lessons are then exemplified in the context of our 
proposed Quantum Irregular Convolutional Code (QIRCC), 
which constitutes the outer component of a concatenated 
quantum code. The proposed QIRCC can be dynamically 
adapted to match any given inner code using EXIT charts, 
hence achieving a performance close to the hashing bound. It 
is demonstrated that our QIRCC-based optimized design is 
capable of operating within 0.4 dB of the noise limit. 

Keywords—Quantum Error Correction, Turbo Codes, EXIT 
Charts, Hashing Bound. 
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BCH 

Bose-Chaudhuri-Hocquenghem 

BIBD 

Balanced Incomplete Block Designs 

BSC 

Binary Symmetric Channel 

CCC 

Classical Convolutional Code 

CSS 

Calderbank-Shor-Steane 

EA 

Entanglement-Assisted 

EXIT 

EXtrinsic Information Transfer 

IRCC 

IRregular Convolutional Code 

LDGM 

Low Density Generator Matrix 

LDPC 

Low Density Parity Check 

MAP 

Maximum A Posteriori 

MI 

Mutual Information 

PCM 

Parity Check Matrix 

QBER 

Qubit Error rate 

QC 

Quasi-Cyclic 

QCC 

Quantum Convolutional Code 

QECC 

Quantum Error Correction Code 

QIRCC 

Quantum IRregular Convolutional Code 

QLDPC 

Quantum Low Density Parity Check 

QSC 

Quantum Stabilizer Code 

QTC 

Quantum Turbo Code 

RX 

Receiver 
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SISO 

Soft-In Soft-Out 

SNR 

Signal-to-Noise Ratio 

TX 

Transmitter 

WER 

Word Error Rate 


I. Introduction 

The laws of quantum mechanics provide a promising solu¬ 
tion to our quest for miniaturization and increased processing 
power, as implicitly predicted by Moore’s law formulated 
four decades ago hi. This can be attributed to the inherent 
parallelism associated with the quantum bits (qubits). More 
explicitly, in contrast to the classical bits, which can either 
assume a value of 0 or 1, qubits can exist in a superposition 
of the two states^]. Consequently, while an A r -bit classical 
register can store only a single value, an TV -qubit quantum 
register can store all the 2 N states concurrently allowing 
parallel evaluations of certain functions with regular global 
structure at a complexity cost that is equivalent to a single 
classical evaluation a, sa, as illustrated in Fig. Q] There¬ 
fore, as exemplified by Shor’s factorization algorithm a and 
Grover’s search algorithm a, quantum-based computation is 
capable of solving certain complex problems at a substantially 
lower complexity, as compared to its classical counterpart. 
From the perspective of telecommunications, this quantum 
domain parallel processing seems to be a plausible solution 
for the massive parallel processing required for achieving joint 
optimization in large-scale communication systems, e.g. quan¬ 
tum assisted multi-user detection ||4), (9j, HOl and quantum- 
assisted routing optimization for self-organizing networks ED. 
Furthermore, quantum-based communication is capable of sup¬ 
porting secure data dissemination, where any ‘measurement’ 
or ‘observation’ by an eavesdropper destroys the quantum 


’The superimposed state of a qubit may be represented as | ip) = a|0) + 
/3|1), where | ) is called Dirac notation or Ket (2), which is a standard notation 
for states in quantum physics, while a and (3 are complex numbers with 
\a\ * 2 + |/3| 2 = 1. More specifically, a qubit exists in a continuum of states 
between |0) and |1) until it is ‘measured’ or ‘observed’. Upon ‘measurement’ 
it collapses to the state |0) with a probability of |a| 2 and |1) with a probability 
of |/3| 2 . 

2 A single qubit is essentially a vector in the 2-dimensional Hilbert space. 
Consequently, an 7V-qubit composite system, which consists of N qubits, has 
a 2^-dimensional Hilbert space, which is the tensor product of the Hilbert 
space of the individual qubits. The resulting iV-qubit state may be generalized 
as: 

a o |00 . .. 0) + ai|00 . .. 1> + • • • + a 2 jv_! |H ... 1> 5 
2^-1 

where a* E C and \ a i\ 2 — 1- 

i=0 
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Fig. 1. Quantum Parallelism: Given a function f(x), which has a regular global structure such that f(x) : {0, l } 2 —»• {0, l} 2 , a classical system requires four 
evaluations to compute f(x) for all possible x E {00, 01,10,11}. By contrast, since a 2-qubit quantum register can be in a supeiposition of all the four states 
concurrently, i.e. \ip) = ao| 00 ) + ai| 01 ) + ol ^| 10 ) + 0 : 3111 ), quantum computing requires only a single classical evaluation to yield the outcome, which is 
also in a superposition of all the four possibilities, i.e. ao|/(00)) + o;i |/(01)) + £* 21 /( 10 )) +q:3 |/(11)). However, it is not possible to read all the four states 
because the quantum register collapses to one of the four superimposed states upon measurement. Nevertheless, we may manipulate the resultant supeiposition 
of the four possible states before observing the quantum register for the sake of determining a desired property of the function, as in ED i 


entanglement hence intimating the parties concerned 0, 
na. Quantum-based communication has given rise to a new 
range of security paradigms, which cannot be created using 
a classical communication system. In this context, quantum 
key distribution techniques M, ED. quantum secure direct 
communication ed, m and the recently proposed uncon¬ 
ditional quantum location verification ED are of particular 
significance. 

Unfortunately, a major impediment to the practical real¬ 
ization of quantum computation as well as communication 
systems is quantum noise, which is conventionally termed 
as ‘decoherence’ (loss of the coherent quantum state). More 
explicitly, decoherence is the undesirable interaction of the 
qubits with the environment fl9l , f20l . It may be viewed as 
the undesirable entanglement of qubits with the environment, 
which perturbs the fragile superposition of states, thus leading 
to the detrimental effects of noise. The overall decoherence 
process may be characterized either by bit-flips or phase- 
flips or in fact possibly both, inflicted on the qubits fl9l , as 
depicted in Fig. m The longer a qubit retains its coherent state 
(this period is known as the coherence time), the better. This 
may be achieved with the aid of Quantum Error Correction 
codes (QECCs), which also rely on the peculiar phenomenon 


3 Two qubits are said to be entangled if they cannot be decomposed 
into the tensor product of the constituent qubits. Let us consider the state 
\ip) = a:|00) + /3|11), where both a and (3 are non-zero. It is not possible to 
decompose it into two individual qubits because we have: 

a|00)+/3|ll> ^ (ai|O)+/3i|l))0(a 2 |O) + /5 2 |l», 

for any choice of oli and subject to normalization. Consequently, a peculiar 
link exists between the two qubits such that measuring one qubit also collapses 
the other, despite their spatial separation. More specifically, if we measure the 
first qubit of we may obtain a | 0 ) with a probability of \a\ 2 and a | 1 ) with 
a probability of |/5| 2 . If the first qubit is found to be |0), then the measurement 
of the second qubit will definitely be |0). Similarly, if the first qubit is 11), 
then the second qubit will also collapse to |1). This mysterious correlation 
between the two qubits, which doesn’t exist in the classical world, is called 
entanglement. It was termed ‘spooky action at a distance’ by Einstein □3 
4 A qubit may be realized in different ways, e.g. two different photon 
polarizations, different alignments of a nuclear spin, two electronic levels of 
an atom or the charge/current/energy of a Josephson junction. 



Fig. 2. Quantum decoherence characterized by bit and phase flips. The 
vertical polarization represents the state | 1 ), while the horizontal polarization 
represents the state | 0 ). 


of entanglement - hence John Preskill eloquently pointed out 
that we are “fighting entanglement with entanglement” HD- 
More explicitly, analogously to the classical channel coding 
techniques, QECCs rectify the impact of quantum noise (bit 
and phase flips) for the sake of ensuring that the qubits retain 
their coherent quantum state with a high fidclityO, thus in effect 
beneficially increasing the coherence time of the unperturbed 
quantum state. This has been experimentally demonstrated 
in E3H35). 

Similar to the family of classical error correction codes 11251 . 
ED, which aim for operating close to Shannon’s capacity 
limit, QECCs are designed to approach the quantum capac¬ 
ity |28H30), or more specifically the hashing bound, which 
is a lower bound of the achievable quantum capacity. A 
significant amount of work has been carried out over the 
last few decades to achieve this objective. However, the field 
of quantum error correction codes is still not as mature as 
that of their classical counterparts. Recently, inspired by the 

5 Fidelity is a measure of closeness of two quantum states E). 
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family of classical near-capacity concatenated codes, which 
rely on iterative decoding schemes, e.g. ED, ED, substantial 
efforts have been invested in If33l - f35l to construct comparable 
quantum codes. In the light of this increasing interest in 
conceiving hashing bound approaching concatenated quantum 
code design principles, the contributions of this paper are: 

1) We survey the evolution towards constructing hash¬ 
ing bound approaching concatenated quantum codes 
with the aid of EXtrinsic Information Transfer (EXIT) 
charts. More specifically, to bridge the gap between 
the classical and quantum channel coding theory, we 
provide insights into the transformation from the family 
of classical codes to the class of quantum codes. 

2) We propose a generically applicable structure for 
Quantum Irregular Convolutional Codes (QIRCCs), 
which can be dynamically adapted to a specific ap¬ 
plication scenario for the sake of facilitating hashing 
bound approaching performance. This is achieved with 
the aid of the EXIT charts ofm 

3) More explicitly, we provide a detailed design example 
by constructing a 10 -subcode QIRCC and use it as an 
outer code in concatenation with the non-catastrophic 
and recursive inner convolutional code of m, m . 
Our QIRCC-based optimized design outperforms both 
the design of m, as well as the exhaustive-search 
based optimized design ofm 

This paper is organized as follows. We commence by 
outlining our design objectives in Section [TT] We then provide 
a comprehensive historical overview of QECCs in Section Hill 
We detail the underlying stabilizer formalism in Section HVlbv 
providing insights into constructing quantum stabilizer codes 
by cross-pollinating their design with the aid of the well- 
known classical codes. We then proceed with the design of 
concatenated quantum codes in Section [V] with a special em¬ 
phasis on their code construction as well as on their decoding 
procedure. In Section [Vi] we will detail the EXIT-chart aided 
code design principles, providing insights into the application 
of EXIT charts for the design of quantum codes. We will then 
present our proposed QIRCC design example in Section I VIII 
followed by our simulation results in Section IVIIII Finally, our 
conclusions and design guidelines are offered in Section HXl 


II. Design Objectives 

Meritorious families of quantum error correction codes can 
be derived from the known classical codes by exploiting 
the underlying quantum-to-classical isomorphism, while also 
taking into account the peculiar laws of quantum mechanics. 
This transition from the classical to the quantum domain must 
address the following challenges 02 : 

• No-Cloning Theorem: Most classical codes are based 
on the transmission of multiple replicas of the same bit, 
e.g. in a simple rate-1/3 repetition code each information 
bit is transmitted thrice. This is not possible in the quan¬ 
tum domain according to the no-cloning theorem 071 . 
which states that an arbitrary unknown quantum state 


cannot be copied/clonecQ. 

• Continuous Nature of Quantum Errors: In contrast to 
the classical errors, which are discrete with bit-flip being 
the only type of error, a qubit may experience both a bit 
error as well as a phase error or in fact both, as depicted 
in Fig. [2] These impairments have a continuous nature 
and the erroneous qubit may lie anywhere on the surface 
of the Bloch spher^D 

• Qubits Collapse upon Measurement: ‘Measurement’ 
of the received bits is a vital step representing a hard- 
decision operation in the field of classical error correc¬ 
tion, but this is not feasible in the quantum domain, since 
qubits collapse to classical bits upon measurement. 

In a nutshell, a classical (n, k) binary code is designed to 
protect discrete-valued message sequences of length k by 
encoding them into one of the 2 k discrete codewords of length 
n. By contrast, since a quantum state of k qubits is specified 
by 2 k continuous-valued complex coefficients, quantum error 
correction aims for encoding a fc-qubit state into an n-qubit 
state, so that all the 2 k complex coefficients can be perfectly 
restored (38]. For example, let k = 2, then the 2-qubit 
information word \rp) is given by: 

| i>) = do 100 ) + ai| 01 ) + a 2 | 10 ) + a 3 |ll). ( 1 ) 

Consequently, the error correction algorithm would aim for 
correctly preserving all the four coefficients, i.e. ao, ai, 0:2 and 
03 . It is interesting to note here that although the coefficients 
ao, a 1 , 09 and a 3 are continuous in nature, yet the entire 
continuum of errors can be corrected, if we can correct a 
discrete set of errors, i.e. bit (Pauli-XfJ phase (Pauli-Z) 
as well as both (Pauli-Y) errors inflicted on either or both 
qubits Ifl3l. This is because measurement results in collapsing 


6 No-cloning theorem is a direct consequence of the linearity of transforma- 
tions. Let us assume that U is a copying operation, which maps the arbitrary 
states \i/j) and \(j>) as follows: 

U\fi) = \rp) U\fi) = |0) <g> |0). 

Since the transformation U must be linear, we should have: 

u ( 10) + 10)) = u \ i >) + u \< t >) = 10) ® 10) +10) ® 10). 

However, 

IVO ® IVO + IVO <s> 1 4 >) / (IVO + I VO) ® (IVO + IVO)- 


7 A qubit |VO = a 10) + /5| 1), whose orthogonal basis are |0) and |1), can 
be visualized in 3D as a unique point on the surface of a sphere (with unit 
radius) called Bloch sphere HE 

8 A qubit |VO = a|0)+/3|l) may be represented as in vector notation. 

Consequently, I, X, Y and Z Pauli operators (or gates), which act on a single 
qubit, are defined as follows: 


I 


1 0 
0 1 






where the X, Y and Z operators anti-commute with each other. The output 
of a Pauli operator may be computed using matrix multiplication, e.g.: 

X(aj0) + /3|1>) = (J J) x Q = (f) = j9|0) + ojl). 
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the entire continuum of errors to a discrete set. More explicitly, 
for \ip) of Eq. 0, the discrete error set is as follows: 

{IX, IZ, IY, XI, XX, XZ, XY, ZI, ZX, ZZ, 

ZY, YI, YX, YZ, YY}. (2) 

However, the errors X, Y and Z may occur with varying 
frequencies. In this paper, we will focus on the specific design 
of codes conceived for mitigating the deleterious effects of 
the quantum depolarizing channel, which has been extensively 
investigated in the context of QECCs If38l - ll40l . Briefly, a 
depolarizing channel, which is characterized by the probability 
p, inflicts an error V G Q„ on n qubitfl where each qubit may 
independently experience either a bit flip (X), a phase flip (Z) 
or both (Y) with a probability of p/3. 

An ideal code C designed for a depolarizing channel may be 
characterized in terms of the channel’s depolarizing probability 
p and its coding rate Rq. Here the coding rate Rq is measured 
in terms of the number of qubits transmitted per channel use, 
i.e. we have Rq = k/n, where k and n are the lengths of the 
information word and codeword, respectively. Analogously to 
Shannon’s classical capacity, the relationship between p and 
Rq for the depolarizing channel is defined by the hashing 
bound, which sets a lower limit on the achievable quantum 
canacit\PT. The hashing bound is given by Il34l . P3l : 

Cq(p) = 1 - H 2 (p) - plog 2 (3), (3) 

where H 2 (p) is the binary entropy function. More explicitly, 
for a given p, if a random code C of a sufficiently long 
codeword-length is chosen such that its coding rate obeys 
Rq < Cq(p), then C may yield an infinitesimally low Qubit 
Error Rate (QBER) for a depolarizing probability of p. It must 
be noted here that intuitively a low QBER corresponds to a 
high fidelity between the transmitted and the decoded quantum 
state. More explicitly, for a given value of p, Cq (p) gives the 
hashing limit on the coding rate. Alternatively, for a given 
coding rate Rq , where we have Rq = Cq(p*), p* gives 
the hashing limit on the channel’s depolarizing probability. In 
duality to the classical domain, this may also be referred to 
as the noise limit. An ideal quantum code should be capable 
of ensuring reliable transmission close to the noise limit p*. 
Furthermore, for any arbitrary depolarizing probability p, its 
discrepancy with respect to the noise limit p* may be computed 


9 A single qubit Pauli group Qi is a group formed by the Pauli matrices I. 
X, Y and Z, which is closed under multiplication. Therefore, it consists of 
all the Pauli matrices together with the multiplicative factors ±1 and ±i, i.e. 
we have: 


G i = {±1, ±il, ±X, ±iX, ±Y. ±iY, ±Z, ±iZ}. 

The general Pauli group Q n is an n-fold tensor product of Q\. 

ll1 Quantum codes are inherently degenerate in nature because different 
errors may have the same impact on the quantum state. For example, let 
| ip) = 100) + 111). Both errors IZ and ZI acting on \ip) yield the same 
corrupted state, i.e. (100) — 111)), and are therefore classified as degenerate 
errors. Due to this degenerate nature of the channel errors, the ultimate capacity 
of quantum channel can be higher than that defined by the hashing bound ED, 
E3- However, none of the codes known to date outperform the hashing bound 
at practically feasible frame lengths. 


in decibels (dB) as follows (34): 


Distance from hashing bound = 10 x log 10 



(4) 


Consequently, our quantum code design objective is to min¬ 
imize the discrepancy with respect to the hashing bound, 
thereby yielding a hashing bound approaching code design. 

It is pertinent to mention here the Entanglement-Assisted 
(EA) regime of B4H47), where the entanglement-assisted 
code C is characterized by an additional parameter c. Here 
c is the number of entangled qubits pre-shared between the 
transmitter and the receiver, thus leading to the terminology 
of being entanglement-assistecf/]. It is assumed furthermore 
that these pre-shared entangled qubits are transmitted over a 
noiseless quantum channel. The resultant EA hashing bound 
is given by lf34l . Il48ll : 


Cq(p) = 1 - H 2 (j>) - plog 2 (3) + E, (5) 


where the so-called entanglement consumption rate is E = 
-. Furthermore, the value of c may be varied from 0 to a 
maximum of (n — k). For the family of maximally entangled 
codes associated with c = (n — k ), the EA hashing bound of 
Eq. 0 is reduced to 041 , Il48l : 


Cq(p) = i- 


( 6 ) 


Therefore, the resultant hashing region of the EA communica¬ 
tion is bounded by Eq. 0 and Eq. 0, which is also illustrated 
in Fig. [3] To elaborate a little further, let us assume that the 
desired coding rate is Rq = 0.4. Then, as gleaned from Fig0 
the noise limit for the ‘unassisted’ quantum code is around 
p* = 0.095, which increases to around p* = 0.25 with the aid 
of maximum entanglement, i.e. we have E = 1 — Rq = 0.6. 
Furthermore, 0 < E < 0.6 will result in bearing noise limits 
in the range of 0.095 < p* < 0.25. Let us assume furthermore 
that we design a maximally entangled code C for Rq = 0.4, 
so that it ensures reliable transmission for p < 0.15. Based on 
Eq. 0, the performance of this code (marked with a circle in 
Fig. 0 is around [10 x logl0(^||)] = 2 dB away from the 
noise limit. We may approach the noise limit more closely by 
optimizing a range of conflicting design challenges, which are 
illustrated in the stylized representation of Fig. 0 For example, 
we may achieve a lower QBER by increasing the code length. 
However, this in turn incurs longer delays. Alternatively, we 
may resort to more complex code designs for reducing the 
QBER, which may also be reduced by employing codes having 
lower coding rates or higher entanglement consumption rates, 
thus requiring more transmitted qubits or entangled qubits. 
Hence striking an appropriate compromise, which meets these 
conflicting design challenges, is required. 


11 A quantum code without pre-shared entanglement, i.e. c = 0, may be 
termed as an unassisted quantum code. EA quantum codes will be discussed 
in detail in Section Hv-El 
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Fig. 3. Unassisted and EA hashing bounds characterized by Eq. 
and Eq. {b}, respectively, for the quantum depolarizing channel. The region 
enclosed by these two bounds, which is labeled the hashing region, defines the 
capacity for varying number of pre-shared entangled qubits (c). At Rq =0.4, 
the unassisted hashing bound gives a noise limit of p* = 0.095, while the 
maximally entangled hashing bound increases the limit to p* = 0.25. The 
circle represents a maximally entangled code with Rq =0.4, which ensures 
reliable transmission for p < 0.15, thus operating at a distance of 2 dB from 
the noise limit. 

Coding Rate 



Fig. 4. Stylized illustration of the conflicting design challenges, which are 
involved in the design of quantum codes. 


III. Historical Overview of Quantum Error 
Correction Codes 

A major breakthrough in the field of quantum information 
processing was marked by Shor’s pioneering work on quan¬ 
tum error correction codes, which dispelled the notion that 
conceiving QECCs was infeasible due to the existence of the 
no-cloning theorem. Inspired by the classical 3-bit repetition 
codes, Shor conceived the first quantum code in his seminal 
paper OH, which was published in 1995. The proposed code 
had a coding rate of 1/9 and was capable of correcting only 
single qubit errors. This was followed by Calderbank-Shor- 
Steane (CSS) codes, invented independently by Calderbank 
and Shor 149) as well as by Steane ED, ED," which facilitated 


the design of good quantum codes from the known classical 
binary linear codes. More explicitly, CSS codes may be defined 
as follows: 

An [n, k\ — fca] CSS code, which is capable of correcting 
t bit errors as well as phase errors, can be constructed from 
classical linear block codes C\(n , k\) and C 2 (n , kf), if C '2 C 
Ci and both C\ as well as the c/maQ of C 2 , i.e. C 2 , can 
correct t errors. Here, C\ is used for correcting bit errors, 
while C '2 is used for phase-error correction. 

Therefore, with the aid of CSS construction, the overall 
problem of finding good quantum codes was reduced to finding 
good dual-containin or self-orthogonal classical codes. Fol¬ 
lowing these principles, the classical [7,4, 3] Hamming code 
was used to design a 7-qubit Steane code OT l having a coding 
rate of 1/7, which is capable of correcting single isolated errors 
inflicted on the transmitted codewords. Finally, Laflamme et 
al. l52l and Bennett et al. K31 independently proposed the 
optimal single error correcting code in 1996, which required 
only 4 redundant qubits. 

Following these developments, Gottesman formalized the 
notion of constructing quantum codes from the classical binary 
and quaternary codes by establishing the theory of Quantum 
Stabilizer Codes (QSCs) lf53l in his Ph.D thesis 154) . In 
contrast to the CSS construction, the stabilizer formalism 
defines a more general class of quantum codes, which imposes 
a more relaxed constraint than the CSS codes. Explicitly, 
the resultant quantum code structure can either assume a 
CSS or a non-CSS (also called unrestricted) structure, but it 
has to meet the symplectic product criteriorQ More specif¬ 
ically, stabilizer codes constitute a broad class of quantum 
codes, which subsumes CSS codes as a subclass and has 
undoubtedly provided a firm foundation for a wide variety 
of quantum codes developed, including for example quantum 
Bose-Chaudhuri-Hocquenghem (BCH) codes l55l - l58) . quan¬ 
tum Reed-Solomon codes (59), (60) , Quantum Low Density 
Parity Check (QLDPC) codes 138) . I6D - I163) . Quantum Con¬ 
volutional Codes (QCCs) |64) - Il67) . Quantum Turbo Codes 
(QTCs) 133) . ||39l as well as quantum polar codes KOI . 
(68) . 169) . These major milestones achieved in the history of 
quantum error correction codes are chronologically arranged 
in Fig. [3 Let us now look deeper into the development of 
QCCs, QLDPC codes and QTCs, which have been the prime 
focus of most recent research both in the classical as well as 
in the quantum domain. 

The inception of QCCs dates back to 1998. Inspired by the 
higher coding efficiencies of Classical Convolutional Codes 
(CCCs) as compared to the comparable block codes and the 
low latency associated with the online encoding and decoding 
of CCCs 1701 , Chau conceived the first QCC in 171) . He also 
generalized the classical Viterbi decoding algorithm for the 
class of quantum codes in 172). but he overlooked some crucial 


12 If G and H are the generator and parity check matrices for any linear 
block code C, then its dual code C 1 - is a unique code with H T and G T as 
the generator and parity check matrices respectively. 

13 Code C with parity check matrix H is said to be dual-containing if it 
contains its dual code C - 1 -, i.e. C 1 - C C and HH T = 0. 

14 Further details are given in Section liv-cl 
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Author(s) 

Coding Efficiency 

Decoding Complexity 

Ollivier and Tillich [64]], 

65] 

Low 

Moderate 

Almeida and Palazzo [73 


Moderate 

Moderate 

Forney et al. [ 66 ], |67| 

High 

Low 


TABLE I. Comparison of the Quantum Convolutional Code (QCC) structures. 


Year 

Author(s) 

Contribution 

1998 

Chau l7ll 

The first QCCs were developed. Unfortunately, some important encoding/decoding aspects were ignored. 

1999 

Chau |72| 

Classical Viterbi decoding algorithm was generalized to the quantum domain. However, similar to |71[, some crucial 
encoding/decoding aspects were overlooked. 

2003 

Ollivier and Tillich 1641, 1651 

Stabilizer-based convolutional codes and their maximum likelihood decoding using the Viterbi algorithm were revisited 
to overcome the deficiencies of [711, [721. Failed to provide better performance or decoding complexity than the 
comparable block codes. 

2004 

Almeida and Palazzo 1731 

Shor-type concatenated QCC was conceived and classical syndrome trellis was invoked for decoding. A high coding 

efficiency was achieved at the cost of a relatively high encoding complexity. 

2005 

Forney et al. [661, [671 

Unrestricted and CSS-type QCCs were derived from arbitrary classical self-orthogonal F 4 and F 2 CCCs, respectively, 
yielding a higher coding efficiency as well as a lower decoding complexity than the comparable block codes. 

2005 

Grassl and Rotteler [74], [751 

Conceived a new construction for QCCs from the classical self-orthogonal product codes. 

2007 

Aly et al. [76 


Algebraic QCCs dervied from BCH codes. 

2008 

Aly et al. [77 


Algebraic QCCs constructed from Reed-Solomon and Reed-Muller Codes. 

2013 

Pelchat and Poulin [781 

Degenerate Viterbi decoding was conceived, which runs the MAP algorithm over the equivalent classes of degenerate 
errors, thereby improving the performance. 


TABLE II. Major contributions to the development of Quantum Convolutional Codes (QCCs). 


Stabilizer codes 1531^1541 
Quantum Reed-Solomon codes H60I 



Fig. 5. Major milestones achieved 
codes. 


Shor’s code [13 

CSS codes ED-123, 5-qubit code t43l . 1521 
Quantum BCH codes ESI, E3i 

Quantum Reed-Muller codes ED 

Quantum LDPC codes ED 

Quantum convolutional codes II641 

Quantum turbo codes m 
Quantum polar codes 14011 
i the history of quantum error correction 


encoding and decoding aspects. Later, Ollivier et al. {64), 
|[65l revisited the class of stabilizer-based convolutional codes. 
Similar to the classical Viterbi decoding philosophy, they also 
conceived a look-up table based quantum Viterbi algorithm 
for the maximum likelihood decoding of QCCs, whose com¬ 
plexity increases linearly with the number of encoded qubits. 
Ollivier et al. also derived the corresponding online encoding 
and decoding circuits having complexity which increased 
linearly with the number of encoded qubits. Unfortunately, 
their proposed rate-1/5 single-error correcting QCC did not 
provide any performance or decoding complexity gain over the 
rate-1/5 single-error correcting block code of (52). Pursuing 
this line of research, Almeida et al. (73) constructed a rate- 
1 /4 single-error correcting Shor-type concatenated QCC from 
a classical CC(2,1,2) and invoked the classical syndrome- 
based trellis decoding for the quantum domain. Hence, the 
proposed QCC had a higher coding rate than the QCC of (64) . 
(65). However, this coding efficiency was achieved at the cost 
of a relatively high encoding complexity associated with the 


concatenated trellis structure. It must be pointed out here that 
the pair of independent trellises used for decoding the bit and 
phase errors impose a lower complexity than a large joint 
trellis would. Finally, Forney et al. (66), (67) designed rate- 
(n — 2) j'n QCCs comparable to their classical counterparts, 
thus providing higher coding efficiencies than the comparable 
block codes. Forney et al. na, (67) achieved this by invoking 
arbitrary classical self-orthogonal rate-1 /n F 4 -linear and F 2 - 
linear convolutional codes for constructing unrestricted and 
CSS-type QCCs, respectively. Forney et al. (66), (67) also 
conceived a simple decoding algorithm for single-error correct¬ 
ing codes. Both the coding efficiency and the decoding com¬ 
plexity of the aforementioned QCC structures are compared in 
Table J] Furthermore, in the spirit of finding new constructions 
for QCCs, Grassl et al. G3), (75) constructed QCCs using 
the classical self-orthogonal product codes, while Aly et al. 
explored various algebraic constructions in (25) and G3, 
where QCCs were derived from classical BCH codes and 
Reed-Solomon and Reed-Muller codes, respectively. Recently, 
Pelchat and Poulin made a major contribution to the decoding 
of QCCs by proposing degenerate Viterbi decoding ED. which 
runs the Maximum A Posteriori (MAP) algorithm (27) over 
the equivalent classes of degenerate errors, thereby improving 
the attainable performance. The major contributions to the 
development of QCCs are summarized in Table [II] 

Although convolutional codes provide a somewhat better 
performance than the comparable block codes, yet they are 
not powerful enough to yield a capacity approaching perfor¬ 
mance, when used on their own. Consequently, the desire 
to operate close to the achievable capacity of Fig. [3] at an 
affordable decoding complexity further motivated researchers 
to design beneficial quantum counterparts of the classical 
LDPC codes (79) , which achieve information rates close to 
the Shannonian capacity limit with the aid of iterative decoding 
schemes. Furthermore, the sparseness of the LDPC matrix is of 
particular interest in the quantum domain, because it requires 
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Year 

Author(s) 

Code Type 

Contribution 



2001 

Postol B61 j 

Non-dual- 
containing CSS 

The first example of QLDPC code constructed from a finite geometry based classical 
code. A generalized formalism for constructing QLDPC codes from the corresponding 
classical codes was not developed. 



2004 

Mackay et al. 1381 

Dual-containing 

CSS 

Various code structures, e.g. bicycle codes and unicycle codes, were conceived for 

constructing QLDPC codes from classical dual-containing LDPC codes. Performance 
impairment due to the presence of unavoidable length-4 cycles was first pointed out 
in this work. Minimum distance of the resulting codes was upper bounded by the 
row weight. 



2005 

Lou et al. |85|, |86| 

Non-dual- 
containing CSS 

The generator and PCM ol classical LDGM codes were exploited for constructing 
CSS codes. An increased decoding complexity was imposed and the codes had an 
upper bounded minimum distance. 


Code Construction 

2007 

Mackay 1801 

Dual-containing 

CSS 

Cayley graph-based QLDPC codes were proposed, which had numerous length-4 
cycles. 

QLDPC 



Camara et al. |63| 

Non-CSS 

QLDPC codes derived from classical self-orthogonal quaternary LDPC codes were 
conceived, which failed to outperform MacKay’s bicycle codes. 



Hagiwara et al. 

HD 

Non-dual- 
containing CSS 

Quasi-cyclic QLDPC codes were constructed using a pair of quasi-cyclic LDPC 
codes, which were found using algebraic combinatorics. The resultant codes had 
at least a girth of 6, but they failed to outperform MacKay’s constructions given 
in (S) 



1008 

Aly et al. [831 

Dual-containing 

CSS 

QLDPC codes were constructed from finite geometries, which failed to outperform 

Mackay’s bicycle codes. 




Djordjevic |&jj 

Dual-containing 

CSS 

BIBDs were exploited to design QLDPC codes, which failed to outperform Mackay’s 
bicycle codes. 



2010 

Tan et al. |91| 

Non-CSS 

Several systematic constructions tor non-CSS QLDPC codes were proposed, four 
of which were based on classical binary quasi-cyclic LDPC codes, while one was 
derived from classical binary LDPC-convolutional codes. These code designs failed 
to outperform Mackay’s bicycle codes. 



2011 

Couvreur et al. 

ED. ED 

Dual-containing 

CSS 

Cayley graph-based QLDPC codes of 1801 were further investigated. The lower bound 
on the minimum distance of the resulting QLDPC was logarithmic in the code length, 
but this was achieved at the cost of an increased decoding complexity. 




Kasai B51. T891 

Non-dual- 
containing CSS 

Quasi-cyclic QLDPC codes of 1871 were extended to non-binary constructions, which 
outperformed Mackay’s bicycle codes at the cost of an increased decoding complexity. 
Performance was still not at par with the classical LDPC codes and minimum distance 
was upper bounded. 




Hagiwara et al. 

HD 

Non-dual- 
containing CSS 

Spatially-coupled QC-QLDPC codes were developed, which outperformed the ‘non- 

coupled’ design of 1871 at the cost of a small coding rate loss. Performance was 
similar to that of 1881, 1891, but larger block lengths were required. 


Decoding 

2008 

Poulin et al. |951 


Heuristic methods were developed to alleviate the performance degradation caused 
by unavoidable length-4 cycles and symmetric degeneracy error. 



2012 

Wang et al. 1961 


Feedback mechanism was introduced in the context of the heuristic methods of 1951 
to further improve the performance. 



20011 

Poulin et al. |33|, 

Non-CSS 

QTCs were conceived based on the interleaved serial concatenation of QCCs. QTCs 




ED 


are free from the decoding issue associated with the length-4 cycles and they offer 


Code Construction 




a wider range of code parameters. Degenerate iterative decoding algorithm was also 
proposed. Unfortunately, QTCs have an upper bounded minimum distance. 

QTC 


2014 

Babar et al. |35] 

Non-CSS 

To dispense with the time-consuming Monte Carlo simulations and to facilitate the 





design of hashing bound approaching QTCs, the application of classical non-binary 
EXIT charts of |97| was extended to QTCs. 


Decoding 

2014 

Wilde et al. 1341 


The iterative decoding algorithm of 1331, 1391 failed to yield performance similar 

to the classical turbo codes. The decoding algorithm was improved by iteratively 
exchanging the extrinsic rather than the a posteriori information. 


TABLE III. Major contributions to the development of iterative quantum codes. 


only a small number of interactions per qubit during the error 
correction procedure, thus facilitating fault-tolerant decoding. 
Moreover, this sparse nature also makes QLDPC codes highly 
degenerate. In this context, Postol (61 j conceived the first ex¬ 
ample of a non-dual-containing CSS-based QLDPC code from 
a finite geometry based classical LDPC in 2001. However, 
he did not present a generalized formalism for constructing 
QLDPC codes from the corresponding classical codes. Later, 
Mackay et al. |38l proposed various code structures (e.g. bicy¬ 
cle codes and unicycle codes) for constructing QLDPC codes 
from the family of classical dual-containing LDPC codes. 
Among the proposed constructions, the bicycle codes were 
found to exhibit the best performance. It was observed that 
unlike good classical LDPC codes, which have at most a single 
overlap between the rows of the Parity Check Matrix (PCM), 
dual-containing QLDPC codes must have an even number of 


overlaps. This in turn results in many unavoidable length-4 
cycles, which significantly impair the attainable performance 
of the message passing decoding algorithm. Furthermore, the 
minimum distance of the proposed codes was upper bounded 
by the row weight. Additionally, Mackay et al. also proposed 
the class of Cayley graph-based dual-containing codes in (80), 
which were further investigated by Couvreur et al. in ED. 
(82). Cay ley-graph based constructions yield QLDPC codes 
whose minimum distance has a lower bound, which is a 
logarithmic function of the code length, thus the minimum 
distance can be improved by extending the codeword (or 
block) length, albeit again, only logarithmically. However, this 
is achieved at the cost of an increased decoding complexity 
imposed by the row weight, which also increases logarith¬ 
mically with the code length. Aly et al. contributed to these 
developments by constructing dual-containing QLDPC codes 






































from finite geometries in (83) , while Djordjevic exploited the 
Balanced Incomplete Block Designs (BIBDs) in (84) . albeit 
neither of these provided any gain over Mackay’s bicycle 
codes. Furthermore, Lou et al. (85) . 11861 invoked the non- 
dual-containing CSS structure by using both the generator 
and the PCM of classical Low Density Generator Matrix 
(LDGM) based codes. Unfortunately, the proposed LDGM 
based constructions also suffered from length-4 cycles, which 
in turn required a modified Tanner graph and code doping for 
decoding, thereby imposing a higher decoding complexity. The 
only exceptions to length-4 cycles were constituted by the class 
of Quasi-Cyclic (QC) QLDPC codes conceived by Hagiwara et 
al. Ell, whereby the constituent PCMs of non-dual-containing 
CSS-type QLDPCs were constructed from a pair of QC-LDPC 
codes found using algebraic combinatorics. The resultant codes 
had at minimum girth of 6, but they did not outperform 
MacKay’s bicycle codes conceived in (38) . Hagiwara’s design 
of 187) was extended to non-binary QLDPC codes in (88) . 
(89), which operate closer to the hashing limit than MacKay’s 
bicycle codes. However, having an upper bounded minimum 
distance remains a deficiency of this construction and the non¬ 
binary nature of the code imposes a potentially high decoding 
complexity. Furthermore, the performance was still not at par 
with that of the classical LDPC codes. The concept of QC- 
QLDPC codes was further extended to the class of spatially- 
coupled QC codes in (90), which outperformed the 'noil- 
coup led’ design of (87) at the cost of a small coding rate loss. 
The spatially-coupled QC-QLDPC was capable of achieving a 
performance similar to that of the non-binary QC-LDPC code 
only when its block length was considerably higher. While 
all the aforementioned QLDPC constructions were CSS-based, 
Camara et al. (63) were the first authors to conceive non- 
CSS QLDPC codes. They invoked group theory for deriving 
QLDPC codes from the classical self-orthogonal quaternary 
LDPC codes. Later, Tan et al. (91) proposed several systematic 
constructions for non-CSS QLDPC codes, four of which 
were based on classical binary QC-LDPC codes, while one 
was derived from classical binary LDPC-convolutional codes. 
Unfortunately, the non-CSS constructions of 1551 . I9T1 failed 
to outperform Mackay’s bicycle codes. Since most of the 
above-listed QLDPC constructions exhibit an upper bounded 
minimum distance, topological QLDPC0 were derived from 
Kitaev’s construction in l92l - l94| . Amidst these activities, 
which focused on the construction of QLDPC codes, Poulin et 
al. were the first scientists to address the decoding issues of 
QLDPC codes (95| . As mentioned above, most of the QLDPC 
codes consist of unavoidable length-4 cycles. In fact, when 
QLDPC codes are viewed in the quaternary formalism, i.e. 
GF(4), then they must have length-4 cycles, which emerge 
from the symplectic product criterion. These short cycles 
erode the performance of the classic message passing decoding 
algorithm. Furthermore, the classic message passing algorithm 
does not take into account the degenerate nature of quantum 
codes, rather it suffers from it. This is known as the ‘symmetric 
degeneracy error’. Hence, Poulin et al. proposed heuristic 
methods in (95) to alleviate the undesired affects of having 


'’Topological code structures are beyond the scope of this paper. 


short cycles and symmetric degeneracy error, which were 
further improved in (96). The major contributions made in 
the context of QLDPC codes are summarized in Table [HU 
while the most pr omisin g QLDPC construction methods are 
compared in Table II vfFl 

Pursuing further the direction of iterative code structures, 
Poulin et al. conceived QTCs in (33) . (39), based on the 
interleaved serial concatenation of QCCs. Unlike QLDPC 
codes, QTCs offer a complete freedom in choosing the code 
parameters, such as the frame length, coding rate, constraint 
length and interleaver type. Moreover, their decoding is not 
impaired by the presence of length-4 cycles associated with 
the symplectic criterion. Furthermore, in contrast to QLDPC 
codes, the iterative decoding invoked for QTCs takes into 
account the inherent degeneracy associated with quantum 
codes. However, it was found in (33) . (39), (98) that the 
constituent QCCs cannot be simultaneously both recursive and 
noncatastrophic. Since the recursive nature of the inner code 
is essential for ensuring an unbounded minimum distanced, 
whereas the noncatastrophic nature is a necessary condition to 
be satisfied for achieving decoding convergence to a vanish¬ 
ingly low error rate, the QTCs designed in (33), (39l had a 
bounded minimum distance. The QBER performance curves 
of the QTCs conceived in [33], (39) also failed to match 
the classical turbo codes. This issue was dealt with in ea, 
where the quantum turbo decoding algorithm of (33) was 
improved by iteratively exchanging the extrinsic rather than 
the a posteriori information. Furthermore, in (33) . (34) . (39l . 
the optimal components of QTCs were found by analyzing 
their distance spectra, followed by extensive Monte Carlo 
simulations for the sake of determining the convergence thresh¬ 
old of the resultant QTC. In order to circumvent this time- 
consuming approach and to facilitate the design of hashing 
bound approaching QTCs, the application of classical non¬ 
binary EXIT charts (97) was extended to QTCs in (351 . An 
EXIT-chart aided exhaustive-search based optimized QTC was 
also presented in (35) . The major contributions made in the 
domain of quantum turbo codes are summarized in Table Hill 

Some of the well-known classical codes cannot be imported 
into the quantum domain by invoking the aforementioned 
stabilizer-based code constructions because the stabilizer codes 
have to satisfy the stringent symplectic product criterion. 
This limitation was overcome in (44)-(4T| with the notion 
of EA quantum codes, which exploit pre-shared entanglement 
between the transmitter and receiver. Later, this concept was 
extended to numerous other code structures, e.g. EA-QLDPC 
code (99), EA-QCC (TOO), EA-QTC (34), (36) and EA-polar 
codes 111011 1. In [34|, (36), it was also found that entanglement- 
assisted convolutional codes may be simultaneously both re¬ 
cursive as well as non-catastrophic. Therefore, the issue of 
bounded minimum distance of QTCs was resolved with the 
notion of entanglement. Furthermore, EA-QLDPC codes are 
free from length-4 cycles in the binary formalism, which in 


16 The second column indicates ‘short cycles’ in the binary formalism. Recall 
that all QLDPC codes must have short cycles in the quaternary formalism. 

17 Unbounded minimum distance of a code implies that its minimum 
distance increases almost linearly with the interleaver length. 
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Code Construction 

Short Cycles 

Minimum Distance 

Delay 

Decoding Complexity 

Bicycle codes |38|| 

Yes 

Upper Bounded 

Standard 

Standard 

Cayley-graph based codes [80J—|82J 

Yes 

Increases with the code length 

Standard 

Increases with the code length 

LDGM-based codes |85|, 186| 

Yes 

Upper Bounded 

Standard 

High 

Non-binary quasi-cyclic codes [88jj, [89| 

No 

Upper Bounded 

Standard 

High 

Spatially-coupled quasi-cyclic codes 19011 

No 

Upper Bounded 

High 

High 


TABLE IV. Comparison of the Quantum Low Density Parity Check (QLDPC) code structures. 


- First EA-QECC constructed PHI 


- EA stabilizer formalism El-El 


EA quantum LDPC codes 
EA quantum convolutional codes |QT 
EA quantum turbo codes on, on 
EA polar codes mm 


Fig. 6. Major milestones achieved in the history of entanglement-assisted 
quantum error correction codes. 


turn results in an impressive performance similar to that of the 
corresponding classical LDPC codes. Hence, the concept of the 
entanglement-assisted regime resulted in a major breakthrough 
in terms of constructing quantum codes, whose behaviour is 
similar to that of the corresponding classical codes. The major 
milestones achieved in the history of entanglement-assisted 
quantum error correction codes are chronologically arranged 
in Fig. |6] 

In this contribution, we design a novel QIRCC, which may 
be used as an outer component in a QTC, or in fact any 
arbitrary concatenated quantum code structure. Explicitly, the 
proposed QIRCC may be invoked in conjunction with any 
arbitrary inner code (both unassisted as well as entanglement- 
assisted) for the sake of attaining a hashing bound approach¬ 
ing performance with the aid of the EXIT charts of SB- 
More specifically, we construct a 10-subcode QIRCC and 
use it as the outer code in concatenation with the non- 
catastropliic and recursive inner convolutional code of m- 
In contrast to the concatenated code of D34V , which exhibited a 
performance within 0.9 dB of the hashing bound, our QIRCC- 
based optimized design operates within 0.4 dB of the noise 
limit. Furthermore, at a Word Error Rate (WER) of 10 -3 , our 
design outperforms the benchmark designed in m by about 
0.5 dB. Our proposed design also yields a lower error rate 
than the exhaustive-search based optimized design of 135)1 . 


the stabilizer formalism for the sake of ensuring a smooth 
transition from the classical to the quantum domain. 

A. Classical Linear Block Codes 

The stabilizer formalism derives its existence from the 
theory of classical linear block codes. A classical linear block 
code C(n, k) maps fc-bit information blocks onto n-bit code¬ 
words. For small values of k and n, this can be readily achieved 
using a look-up table, which maps the input information blocks 
onto the encoded message blocks. However, for large values 
of k and n, the process may be simplified using an k x n 
generator matrix G as follows: 

x = xG , (7) 

where x and x are row vectors for information and encoded 
messages, respectively. Furthermore, G may be decomposed 
as: 

G=(I k \P), (8) 

where Ik is a (A: x fc)-element identity matrix and P is a k x 
(n — fc)-element matrix. This in turn implies that the first k 
bits of the encoded message are information bits, followed by 
{n — k) parity bits. 

At the decoder, syndrome-based decoding is invoked, which 
determines the position of the channel-induced error using the 
observed syndromes rather than directly acting on the received 
codewords. More precisely, each generator matrix is associated 
with an (n — k) x n-element PCM H which is given by: 

H = ( P T \I n -k) , (9) 

and is defined such that x is a valid codeword only if, 

xH t = 0. (10) 

For a received vector y = x + e, where e is the error incurred 
during transmission, the error syndrome of length (n — k) is 
computed as: 

s = yH T = (x + e)H T = xH T + eH T = eH T , (11) 


IV. Stabilizer Formalism 

Most of the quantum codes developed to date owe their 
existence to the theory of stabilizer codes, which allows us 
to import any arbitrary classical binary as well as quaternary 
code to the quantum domain. Unfortunately, this is achieved 
at the cost of imposing restrictions on the code structure, 
which may adversely impact the performance of the code, 
e.g. as in QLDPC codes and QTCs, which was discussed 
in Section m In this section, we will delve deeper into 


which is then used for identifying the erroneous parity bit. 

Let us consider a simple 3-bit repetition code, which makes 
three copies of the intended information bit. More precisely, 
k = 1 and n = 3. It is specified by the following generator 
matrix: 

G=( 1 1 1), (12) 

which yields two possible codewords [111] and [000]. At the 
receiver, a decision may be made on the basis of the majority 
voting. For example, if y = [011] is received, then we may 
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Syndrome (s) 

Index of Error 

p] 

1 

[10] 

2 

[01] 

3 


TABLE V. Single-bit errors along with the corresponding 
SYNDROMES FOR THE PCM OF EQ. ( fTsl . 


conclude that the transmitted bit was 1. Alternatively, we 
may invoke the PCM-based syndrome decoding. According 
to Eq. the corresponding PCM is given by: 

»=(;;;). <i3> 

It can be worked out that yH T = 0 only for the two valid 
codewords [111] and [000]. For all other received codewords, 
at least one of the two syndrome elements is set to 1 , e.g. 
when the first bit is corrupted, i.e. y = [Oil] or [100], s = [11], 
Table [V] enlists all the 1 -bit errors, which may be identified 
using this syndrome decoding procedure. 

This process of error correction using generator and parity 
check matrices is usually preferred due to its compact nature. 
Generally, C(n, k) code, which encodes a A--hit information 
message into an n-bit codeword, would require 2 k n-bit 
codewords. Thus, it would required a total of n2 k bits to 
completely specify the code space. By contrast, the afore¬ 
mentioned approach only requires kn bits of the generator 
matrix. Hence, memory resources are saved exponentially and 
encoding and decoding operations are efficiently implemented. 
These attractive features of classical block linear codes and the 
associated PCM-based syndrome decoding 111021 have led to 
the development of quantum stabilizer codes. 


B. Quantum Stabilizer Codes (QSCs) 

Let us recall from Section [n]that qubits collapse to classical 
bits upon measurement 03. This prevents us from directly 
applying the classical error correction techniques for reliable 
quantum transmission. Inspired by the PCM-based syndrome 
decoding of classical codes, Gottesman J53), (Ml introduced 
the notion of stabilizer formalism, which facilitates the design 
of quantum codes from the classical ones. Analogous to Shor’s 
pioneering 9-qubit code ED, stabilizer formalism overcomes 
the measurement issue by observing the error syndromes with¬ 
out reading the actual quantum information. More specifically, 
QSCs invoke the PCM-based syndrome decoding approach of 
classical linear block codes for estimating the errors incurred 
during transmission. 

Fig. H] shows the general schematic of a quantum com¬ 
munication system relying on a quantum stabilizer code for 
reliable transmission. An [n, /.;] QSC encodes the information 
qubits \ip) into the coded sequence | ip) with the aid of (n — k) 
auxiliary (also called ancilla) qubits, which are initialized to 
the state |0). The noisy sequence \ip) = V\ip), where V is 
the n-qubit channel error, is received at the receiver (RX), 
which engages in a 3-step process for the sake of recovering 
the intended transmitted information. More explicitly, RX 


computes the syndrome of the received sequence \ip) and uses 
it to estimate the channel error V. The recovery operator 1Z 
then uses the estimated error V to restore the transmitted coded 
stream. Finally, the decoder, or more specifically the inverse 
encoder, processes the recovered coded sequence \ip), yielding 
the estimated transmitted information qubits | ip). 

An [n, k] quantum stabilizer code, constructed over a code 
space C, which maps the information word (logical qubits) 
\ip) £ C 2 onto the codeword (physical qubits) \ip) £ C 2 , 
where <C d denotes the d-dimensional Hilbert space, is defined 
by a set of (n — k) independent commuting n-tuple Pauli 
operators gi, for 1 < i < ( n—k ). The corresponding stabilizer 
group Tl contains both gi and all the products of gi for 
1 < i < (n — k) and forms an abelian subgroup of Q n . A 
unique feature of these operators is that they do not change 
the state of valid codewords, while yielding an eigenvalue of 
— 1 for corrupted states. 

Let us now elaborate on this definition of the stabilizer 
code by considering a simple 3-qubit bit-flip repetition code, 
which is capable of correcting single-qubit bit-flip errors. Since 
the laws of quantum mechanics do not permit cloning of the 
information qubit, we cannot encode \ip) to (ip ® t/> 0 ip). 
Instead, the 3-qubit bit-flip repetition code entangles two 
auxiliary qubits with the information qubit such that the basis 
states | 0 ) and | 1 ) are copied thrice in the superposition of 
basis states of the resulting 3-qubit codeword, i.e. |0) and |1) 
are mapped as follows: 

| 0 ) | 000 >, 

| 1 >-H 111 >. (14) 

Consequently, the information word |^>) = a|0) + /3|1) is 
encoded as: 


a|0) + /3|1) —> ct|000) + /3|111). (15) 


The resultant codeword is stabilized by the operators gi = ZZI 
and c /2 = ZIZ. Here the term ‘stabilize’ implies that the valid 
codewords are not affected by the generators g\ and g -2 and 
yield an eigenvalue of + 1 , as shown below: 


9i [M] = a|000) + /3|111) = IV*), 

J 2 pj=a |000)+/3|111) = |^). (16) 

On the other hand, if a corrupted state | ip) is received, then 
the stabilizer generators yield an eigenvalue of — 1 , e.g. let 
|^) = 1100) + /3|011) where V = XII, then we have: 


9 1 


—a| 100 ) — /?| 011 ) = —\ip), 


9 2 


= —a|100) — /3|011) = —|i/i). (17) 


More specifically, the eigenvalue is —1 if the n-tuple Pauli 
error V acting on the transmitted codeword | ip) anti-commutes 
with the stabilizer gi and it is +1 if V commutes with gi. 
Therefore, we have: 


9i 


IV>) 


W)j_ 9iP = 'Pgi 

9 iP = -Vgi , 


( 18 ) 
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Fig. 7. System Model: Quantum communication system relying on a quantum stabilizer code. 


\i>) 

9 il'0> 

S2 IV’> 

Syndrome (s) 

Index of Error 

a|100) + /3|011) 

-1 

-1 

PT] 

1 

a 010) + /3|101> 

-1 

+1 

[10] 

2 

a|001) + /3|110) 

+ 1 

-1 

[01] 

3 


TABLE VI. Single-qubit bit-flip errors along with the 

CORRESPONDING EIGENVALUES FOR 3-QUBIT BIT-FLIP REPETITION CODE. 


where | ijj) = V\ip). Table IvTl enlists the eigenvalues for all 
possible single-qubit bit-flip errors. The resultant ±1 eigen¬ 
value gives the corresponding error syndrome s, which is 0 
for an eigenvalue of +1 and 1 for an eigenvalue of — 1 , as 
depicted in Table [Vi] 

A 3-qubit phase-flip repetition code may be constructed 
using a similar approach. This is because phase errors in the 
Hadamard basis {|+}, |—}} are similar to the bit errors in the 
computational basis {|0), |1)}. More explicitly, the states |+) 
and |—) are defined as: 


1+) — H|0) = 
|-) = H|l) = 


| 0 ) + | 1 ) 
V2 

|Q)-|i) 

V2 


(19) 


where H is a single-qubit Hadamard gate, 
by m- 


H = 


- 1 1 
V2 V 1 


which is given 
( 20 ) 


Therefore, Pauli-Z acting on the states |+) and |—) yields: 


z l+) — l _ )i 

Z|-) = |+), (21) 


which is similar to the operation of Pauli-X on the states |0) 
and 11 ), i.e. we have: 


X|0) = |1), 

X|l) = |0). (22) 

Consequently, analogous to Eq. (114b , a 3-qubit phase-flip 
repetition code encodes | 0 ) and | 1 ) as follows: 


| 0 ) ->■ | + ++}, 

|1>-H-)■ (23) 

Based on Eq. ( l23l >. | ip) is encoded to: 


a|0) + /3|1) —> a\ + ++) + -), (24) 


which is stabilized by the generators gi = XXI and <?2 = 
XIX. Hence, the Hadamard and Pauli-X operators enable a 
quantum code to correct phase errors. This overall transition 
from the classical 3-bit repetition code of Section HV-AI to the 
quantum repetition code is summarized in Fig. [8] 

Furthermore, the stabilizer generators gi constituting the sta¬ 
bilizer group 7~L must exhibit the following two characteristics: 

1) Any two operators in the stabilizer set must com¬ 
mute so that the stabilizer operators can be applied 
simultaneously, i.e. we have: 

9i92\ip) =g2gi\ip)- (25) 

This is because the stabilizer leaves the codeword 
unchanged as encapsulated below: 

9i\^) = I fp)- (26) 

Hence, evaluating the left-hand and right-hand sides 
of Eq. ( l25l > gives: 

gi92\ip) = 5i|^) = IV*), (27) 

and 

92 gi\i>) = 92 ^) = IV’), (28) 

respectively. This further imposes the constraint that the 
stabilizers should have an even number of places with 
different non-identity (i.e. X, Y, or Z) operations. This 
is derived from the fact that the X, Y and Z operations 
anti-commute with one another as shown below: 

XY = iZ, YX = -iZ aXY = YX (29) 

YZ = zX, ZY = -zX —> YZ = -ZY 

ZX = zY, XZ = -zY —»■ ZX = -XZ 


Thus, for example the operators ZZI and XYZ com¬ 
mute, whereas ZZI and YZI anti-commute. 

2) Generators constituting the stabilizer group % are 
closed under multiplication, i.e. multiplication of 
the constituent generators gi yields another generator, 
which is also part of the stabilizer group 'H. For 
example, the full stabilizer group T~L of the 3-qubit bit- 
flip repetition code will also include the operator IZZ, 
which is the product of g\ and gi- 

It must be mentioned here that the Pauli errors which differ 
only by the stabilizer group have the same impact on all the 
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Fig. 8. Transition from the classical to quantum codes. Encoder: In classical codes, the information bit may be copied during the encoding process, e.g. in a 
3-bit repetition code. This is not permissible in the quantum domain (no-cloning theorem). Alternatively, in quantum codes, the information qubit is entangled 
to the auxiliary qubits, e.g. in a 3-qubit bit-flip repetition code. Channel: Only bit errors occur over a classical channel, e.g. a binary symmetric channel with 
the channel crossover probability p. By contrast, qubits may experience a bit or phase error as well as both, e.g. depolarizing channel with a probability p. Since 
phase errors in the Hadamard basis {|+), |—)} are similar to the bit errors in the computational basis {10), 11)}, phase errors may be corrected in the same 
way as the bit errors by exploiting the Hadamard basis. Decoder: In classical codes, the received bits are measured during the decoding process, e.g. in a 3-bit 
repetition code a decision may be made on the basis of majority voting. Unfortunately, qubits collapse upon measurement. Consequently, quantum codes invoke 
the PCM-based syndrome decoding. 


codewords and therefore can be corrected by the same recovery 
operations. This gives quantum codes the intrinsic property 
of degeneracy [78]. More explicitly, the errors V and V' = 
g{P have the same impact on the transmitted codeword and 
therefore can be corrected by the same recovery operation. 
Getting back to our example of the 3-qubit bit-flip repetition 
code, let V = IIX and V' = g{P = ZZX. Both V as well as 
V' corrupt the transmitted codeword of Eq. ([J5| to a|001) + 
ji\ 110). Consequently, V and V' need not be differentiated and 
are therefore classified as degenerate errors. 


C. Pauli-to-Binary Isomorphism 

QSCs may be characterized in terms of an equivalent clas¬ 
sical parity check matrix notation satisfying the commutativity 
constraint of stabilizers ||38l , 111031 given in Eq. ( |25| >. This is 
achieved by mapping the I, X, Y and Z Pauli operators onto 
(F 2 ) 2 as follows: 

1 ^ ( 0 , 0 ), 

X^ (0,1), 

( 1 , 1 ), 

Z->(1,0). (30) 

More explicitly, the (n — k) stabilizers of an [n, k] stabilizer 
code constitute the rows of the binary PCM H, which can be 
represented as a concatenation of a pair of (n — k) x n binary 
matrices H z and H x based on Eq. ( 13Oi l, as given below: 


H = (H Z \H X ). 


Each row of H corresponds to a stabilizer of H, so that the 
/th column of II z and H x corresponds to the /th qubit and a 
binary 1 at these locations represents a Z and X Pauli operator, 
respectively, in the corresponding stabilizer. For the 3-qubit bit- 
flip repetition code, which can only correct bit-flip errors, the 
PCM H is given by: 


( 1 1 

0 

0 

0 

0 

u ° 

1 

0 

0 

0 ) 


(32) 


It must be pointed out here that H z of Eq. ([32} is same as 
the H of the classical repetition code of Eq. dl 3b , yielding the 
same syndrome patterns in Table [V] and Table [VI] 


Let us further elaborate the process by considering the [9,1] 
Shor’s code, which consists of the Pauli-Z as well as the Pauli- 
X operators. The corresponding stabilizer generators are given 
in Table 1 VIII They can be mapped onto the binary matrix H 
as follows: 


H = 


Hi 0 

0 H' x 


(33) 


where we have H z 




/I 1 0 0 0 

0 110 0 
0 0 0 1 1 

0 0 0 0 1 

0 0 0 0 0 

\0 0 0 0 0 


0 0 0 0 \ 

0 0 0 0 
0 0 0 0 
1 0 0 0 ’ 
0 110 
0 0 11 / 


( 34 ) 


(31) 
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91 

92 

93 

94 

95 

96 

97 

98 


Stabilizer 


ZZIIIIIII 

IZZIIIIII 

IIIZZIIII 

IIIIZZIII 

IIIIIIZZI 

IIIIIIIZZ 

XXXXXXIII 

IIIXXXXXX 


TABLE VII. 


HL = 


Stabilizers for 9-qubit Shor’s code. 


1111110 0 0 
0 0 0 1 1 1 1 1 1 


(35) 


With the matrix notation of Eq. (OH . the commutative 
property of stabilizers given in Eq. ({25} is transformed into the 
orthogonality of rows with respect to the symplectic product 
(also called twisted product). If row m is r m = 
where z m and x m are the binary strings for Z and X respec¬ 
tively, then the symplectic product of rows m and m' is given 

by. 



Fig. 9. Family of stabilizer codes. 


Cm * t’m' — (^m * %m' T m' ‘ %m) Ittod 2. (36) 

This symplectic product is zero if there are even number of 
places where the operators (X or Z) in row m and m! are 
different; thus meeting the commutativity requirement. In other 
words, if H is written as H = (H Z \H X ), then the symplectic 
product is satisfied for all the rows only if, 

H z Hj + H x Hj = 0, (37) 

which may be readily verified for the H of Eq. {33}. Con¬ 
sequently, any classical binary codes satisfying Eq. ( 1371 may 
be used to construct QSCs. A special class of these stabilizer 
codes are CSS codes, which are defined as follows; 

An [n,ki — kf\ CSS code, which is capable of correcting t 
bit as well as phase errors, can be constructed from classical 
linear block codes Ci(n,k i) and C 2 (n, /c 2 ), if Ci C C\ and 
both Ci as well as the dual of C 2 , i.e. Q/, can correct t 
errors. 

In CSS construction, the PCM H' z of C\ is used for correct¬ 
ing bit errors, while the PCM H' x of Cf is used for phase-error 
correction. Consequently, the PCM of the resultant CSS code 
takes the form of Eq. ( 1331 . H' z and H' x are now the (n—ki)x n 
and k -2 x n binary matrices, respectively. Furthermore, since 
C '2 C Ci, the symplectic condition of Eq. {37} is reduced 

T 

to H' Z H' X = 0. In this scenario, (n — ki + stabilizers 
are applied to n qubits. Therefore, the resultant quantum 
code encodes ( ki — k 2 ) information qubits into n qubits. 
Furthermore, if Hi = H' the resultant structure is called dual- 

T 

containing (or self-orthogonal) code because H z rH' z = 0, 
which is equivalent to Ci C C\ . Hence, stabilizer codes 
may be sub-divided into various code structures, which are 
summarized in Fig. [9] 

Let us consider the classical (7,4) Hamming code, whose 
PCM is given by: 

(1 1 0 1 1 0 0\ 

H = [l 0 1 1 0 1 0 . (38) 

\o 1 1 1 0 0 1/ 


Since the H of Eq. {38} yields HH T = 0, it is used for 
constructing the dual-containing rate-1/7 Steane code BTl . 

Based on the aforementioned Pauli-to-binary isomorphism, 
a quantum-based Pauli error operator V can be represented 
by the effective classical error pattern P, which is a binary 
vector of length 2 n. More specifically, P is a concatenation of 
n bits for Z errors, followed by another n bits for X errors, 
as depicted in Fig. [lO] An X error imposed on the 1st qubit 
will yield a 0 and a 1 at the 1st and ( n + l)th index of P, 
respectively. Similarly, a Z error imposed on the 1st qubit 
will give a 1 and a 0 at the 1st and (n + l)th index of P, 
respectively, while a Y error on the 1 st qubit will result in a 1 
at both the 1st as well as (n + l)th index of / F3. The resultant 
syndrome is given by the symplectic product of H and P, 
which is equivalent to H(P, : : P Z ) T . Here colon (:) denotes the 
concatenation operation. In other words, the Pauli-X operator 
is used for correcting Z errors, while the Pauli-Z operator 
is used for correcting X errors lfl3l . Thus, the quantum- 
domain syndrome is equivalent to the classical-domain binary 
syndrome and a basic quantum-domain decoding procedure is 
similar to syndrome based decoding of the equivalent classical 
code |38l . However, due to the degenerate nature of quantum 
codes, quantum decoding aims for finding the most likely error 
coset, while the classical syndrome decoding II 10211 finds the 
most likely error. 

Hence, an [n, k] quantum stabilizer code associated with 
(n — k) stabilizers can be effectively modeled using an (n — 
k) x 2 n-elenient classical PCM satisfying Eq. (ITT} . The coding 
rate of the equivalent classical code R c can be determined as 


18 Since a depolarizing channel characterized by the probability p incurs X, 
Y and Z errors with an equal probability of p/3, the effective error-vector 
P reduces to two Binary Symmetric Channels (BSCs), one channel for the Z 
eirors and the other for the X errors. The crossover probability of each BSC 
is given by 2p/3. 
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Fig. 10. Effective classical error P corresponding to the error V imposed 
on an n-qubit frame. 


follows: 


Rc 


2 n — (n — k ) 
2 n 

n + k 


2 n 



^ (1 + Rq) i 


(39) 


where i?Q is its quantum coding rate. Using Eq. (TT9b . the 
coding rate of the classical equivalent of Shor’s rate-1/9 
quantum code is 5/9. 


D. Stabilizer Formalism of Quantum Convolutional Codes 

Quantum convolutional codes are derived from the corre¬ 
sponding classical convolutional codes using stabilizer formal¬ 
ism. This is based on the equivalence between the classical 
convolutional codes and the classical linear block codes with 
semi-infinite length, which is derived below l26l . 

Consider a (2,1, m) classical convolutional code with gen¬ 
erators. 


9 (0) = ( 9 ™, 9 ™..., 9 fi>), 

9 < 1 > = ( 9 S ,, , 9 i 1) ,..., 9 S>)- (40) 


This encoding process can also be represented in matrix nota¬ 
tion by interlacing the generators g >u> and g' 11 and arranging 
them in matrix form as followr-FT. 

/ (0)(1) (0)(1) (0)(1) \ 

90 9 1 • ■ • 9m 

„(0)(1) „(0)(1) „(0)(1) 

„ 90 91 9m 

G = (0)(i) (0)(i) (0)(i) , 

9 0 9 1 • • • 9m 

V ••• ■■ ) 

(44) 

where gf' 1 ^' 1 — ( 9 ^ 9 ^)- The encoding operation of 

Eq. (l42l > is therefore equivalent to. 


v = uG. (45) 

Since the information sequence u is of arbitrary length, G 
is semi-infinite. Furthermore, each row of G is identical to 
the previous row, but is shifted to the right by two places 
(since n = 2). In practice, u has a finite length N. Therefore, 
G has N rows and 2 (m + N) columns for CC(2,1 ,m). For 
CC (n,k,m), G can be generalized as follows: 


G = 


/G 0 G 1 ... G m 
G 0 Gi ... 
Go G, 

V 


G n 


G n 


■/ 


(46) 


where G; is a (k x n) submatrix with entries, 


(9 ( S 

$ ■ 

• sir 1 ’) 


92 } 

92 } ■ 

• sir 1 ’ 

(47) 

WkJ 

9k} ■ 

■ sir") 



The corresponding PCM H can be represented as a semi¬ 
infinite matrix consisting of submatrices Hi with dimensions 
of ( n—k ) x n. For a convolutional code with constraint lengtfFI 
(m + 1), H is given by: 


For an input sequence [it = (uo, tti, M 2 , • ■ •)], the out¬ 
put sequences [i/°( = (vq°\ v^, v^, ■ ■ ■)] and [t/ 1 ) = 
(uq 1 '* , uf ^, v 2 1 "*,... )]are given as follows: 

ri(°) = u © g(°\ 

yGl = u® (41) 


H = 


(Ho 

H\ Ho 
U> //1 Ho 


Hm H m — 1 H m —2 
Hm H m _ 1 


Hm—2 


Ho 


\ 


Ho 


(48) 


where © denotes discrete convolution (modulo 2 ), which 
implies that for alH > 0 we have: 


M 


= U l-i9i = Ul9o + u l-l9l 4-f Ul-m9m\ 


i=0 


(42) 

where j = 0,1 and ui-i = 0 for all l < i. The two encoded 
sequences are multiplexed into a single codeword sequence v 
given by: 


, ( 0 ) ( 1 ) ( 0 ) ( 1 ) ( 0 ) ( 1 ) 
V = (v'o ,v^ J ,vl ,vl ,v k 2 ’,v k 2 


= (■ 


...) 


( 43 ) 


V : 

Therefore, a CCC can be represented as a linear block code 
with semi-infinite block length. Furthermore, if each row of 
the submatrices Hi is considered as a single block and hjj 
is the ith row of the yth block, then H has a block-band 
structure after the first m blocks, whereby the successive 
blocks are time-shifted versions of the first block (j = 0 ) 
and the adjacent blocks have an overlap of m submatrices. 


Milank spaces in the matrix indicate zeros. 

-'Constraint length is the number of memory units (shift registers) plus 1. 
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(m x n + n) 


H-m H-m—1 H m —2 • • • H 0 



H m H m —i H m —2 . . . 

Ho 


-* - m x n - *- 



Fig. 11. Block-band structure of the semi-infinite classical PCM H. 


This has been depicted in Fig. QT| and can be mathematically 
represented as follows: 

hj,i = [CF X ",M, 1 <i<(n-k), 0 < j, (49) 

where 0 jxn is a row-vector with (j x n) zeros. 

As discussed in Section IIV-CI the rows of a classical PCM 
correspond to the stabilizers of a quantum code. Hence, the 
quantum stabilizer group T~L of an [n, k, m] stabilizer convolu¬ 
tional code is given by (65]: 

U = sp{g J}i = I® jn <g) go,*}, 1 < i < (n - k), 0 < j, (50) 


where g :h , is the vth stabilizer of the j th block of the stabilizer 
group 7~L . Furthermore, sp represents a symplectic group, thus 
implying that all the stabilizers g 3 ^ must be independent and 
must commute with each other. 

As proposed by Forney in [ 66 ], (67], CSS-type QCCs can be 
derived from the classical self-orthogonal binary convolution 
codes. Let us consider the rate 1/3 QCC of ( 66 ], (67), which 
is constructed from a binary rate-1/3 CCC with generators: 


G = 



1 

1 

1 

0 

0 

1 

1 

0 

0 

0 

0 

... \ 

(° 

0 

0 

1 

1 

1 

1 

0 

0 

1 

1 

0 



(51) 


In /7-transform notation, these generators are represented as 
(1 + D + D 2 ,l + D 2 , 1). Each generator is orthogonal to 
all other generators under the binary inner product, making 
it a self-orthogonal code. Moreover, the dual C 1 - has the 
capability of correcting 1 bit. Therefore, based on the CSS 
construction, the basic stabilizers of the corresponding single¬ 
error correcting [3,1] QCC are as follows: 


5 0.1 = [XXX, XII, XXI], (52) 

50.2 = [ZZZ,ZII,ZZI], (53) 

Other stabilizers of 'H are the time-shifted versions of these 
basic stabilizers as depicted in Eq. (15 Oi l. 

Let us further consider a non-CSS QCC construction given 
by Forney in ED, El. It is derived from the classical self- 
orthogonal rate-1/3 quaternary (F 4 ) convolutional code C 
having generators (1 + 0,1 + wD, 1 + wD ), where F 4 = 


{0, I.t+tu}. These generators can also be represented as 
follows: 



1 

1 

1 

w 

w 

0 

0 

0 


(° 

0 

0 

1 

1 

1 

1 

w 

W 

-) 


Since all these generators are orthogonal under the Hermitian 
inner product, C is self-orthogonal. Therefore, a [3,1] QCC 
can be derived from this classical code. The basic generators 
5 o,i, for 1 < i < 2, of the corresponding stabilizer group, 
H, are generated by multiplying the generators of Eq. (54} 
with w and w, and mapping 0, w, 1, w onto I, X, Y and Z 
respectively. The resultant basic stabilizers are as follows: 

50.1 = (XXX, XZY), (55) 

50.2 = (ZZZ, ZYX), (56) 

and all other constituent stabilizers of H can be derived using 
Eq. ®. 


E. Entanglement-Assisted Stabilizer Formalism 


Let us recall that the classical binary and quaternary codes 
may be used for constructing stabilizer codes only if they 
satisfy the symplectic criterion of Eq. (37). Consequently, some 
of the well-known classical codes cannot be explored in the 
quantum domain. This limitation can be readily overcome by 
using the entanglement-assisted stabilizer formalism, which 
exploits pre-shared entanglement between the transmitter and 
receiver to embed a set of non-commuting stabilizer generators 
into a larger set of commuting generators. 

Fig. [T2l shows the general schematic of a quantum commu¬ 
nication system, which incorporates an Entanglement-Assisted 
Quantum Stabilizer Code (EA-QSC). An [?r, fc,c] EA-QSC 
encodes the information qubits \tp) into the coded sequence | ip) 
with the aid of (n— k—c) auxiliary qubits, which are initialized 
to the state |0). Furthermore, the transmitter and receiver share 
c entangled qubits (ebits) before actual transmission takes 
place. This may be carried out during the off-peak hours, when 
the channel is under-utilized, thus efficiently distributing the 
transmission requirements in time. More specifically, the state 
|L + ) of an ebit is given by the following Bell state: 


\<t> + ) 


|00}7+flx + \\\)TxRx 

72 


(57) 


where Tx and Rx denotes the transmitter’s and receiver’s 
half of the ebit, respectively. Similar to the superdense coding 
protocol of 11041 . it is assumed that the receiver’s half of the c 
ebits are transmitted over a noiseless quantum channel, while 
the transmitter’s half of the c ebits together with the ( n — k — c ) 
auxiliary qubits are used to encode the intended k information 
qubits into n coded qubits. The resultant ro-qubit codewords 
1+) are transmitted over a noisy quantum channel. The receiver 
then combines his half of the c noiseless ebits with the received 
n-qubit noisy codewords \ip) to compute the syndrome, which 
is used for estimating the error V incurred on the 71 -qubit 
codewords. The rest of the processing at the receiver is the 
same as that in Fig. [7] 
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Fig. 12. System Model: Quantum communication system relying on an entanglement-assisted quantum stabilizer code 


The entangled state of Eq. (l57l > has unique commutativity 
properties, which aid in transforming a set of non-abelian 
generators into an abelian set. The state |</> + ) is stabilized by 
the operators X Tx X Hx and Z Tx Z Rx , which commute with 
each other. Therefore, we havcFI: 

[X. Tx X Rx ,Z Tx Z Rx ] = 0. (58) 

However, local operators acting on either of the qubits anti¬ 
commute, i.e. we have: 

{X Tx ,Z Tx } = {X Rx ,Z Rx } = 0. (59) 


In Eq. (l62l i. the first two generators (i.e. the first and second 
row) anti-commute, while all other generators commute with 
each other. This is because the local operators acting on the 
second qubit in the first two generators anti-commute, while 
the local operators acting on all other qubits in these two 
generators commute. In other words, there is a single index (i.e. 
2) with different non-identity operators. To transform this non- 
abelian set into an abelian set, we may extend the generators of 
Eq. ( l62l > with a single additional qubit, whose local operators 
also anti-commute for the sake of ensuring that the resultant 
extended generators commute. Therefore, we get: 


Therefore, if we have two single qubit operators X Tx and 
Z Tx , which anti-commute with each other, then we can 
resolve the anti-commutativity by entangling another qubit and 
choosing the local operators on this additional qubit such that 
the resultant two-qubit generators (X Tx X Rx and Z Tx Z Rx 
for this case) commute. This additional qubit constitutes the 
receiver half of the ebit. In other words, we entangle an 
additional qubit for the sake of ensuring that the resultant two- 
qubit operators have an even number of places with different 
non-identity operators, which in turn ensures commutativity. 


Let us consider a pair of classical binary codes associated 
with the following PCMs: 


/0 1 0 0\ 

0 0 0 01 
1110’ 
\0 1 1 1 / 


(60) 


and 


H x 


/I 0 1 0\ 
110 1 
10 0 1 ’ 
\1 1 1 0 / 


(61) 


which are used to construct a non-CSS quantum code having 
H = (H Z \H X ). The PCM H does not satisfy the symplectic 
criterion. The resultant non-abelian set of Pauli generators are 
as follows: 


/ X 
X 
Y 

V X 


z 

X 

z 

Y 


X 

I 

Z 

Y 



(62) 


21 [a, b] represents the commutative relation between a and b, while {a, 6} 
denotes the anti-commutative relation. 


/ X Z X I 

XXIX 
Mq - Y Z Z X 

V X Y Y Z 

where the operators to the left of the vertical bar (|) act on the 
transmitted rc-qubit codewords, while those on the right of the 
vertical bar act on the receiver’s half of the ebits. 

V. Concatenated Quantum Codes 

In this section, we will lay out the structure of a concate¬ 
nated quantum code, with a special emphasis on the encoder 
structure and the decoding algorithm. We commence with 
the circuit-based representation of quantum stabilizer codes, 
followed by the system model and then the decoding algorithm. 



A. Circuit-Based Representation of Stabilizer Codes 

Circuit-based representation of quantum codes facilitates the 
design of concatenated code structures. More specifically, for 
decoding concatenated quantum codes it is more convenient 
to exploit the circuit-based representation of the constituent 
codes, rather than the conventional PCM-based syndrome 
decoding. Therefore, in this section, we will review the circuit- 
based representation of quantum codes. This discussion is 
based on 133). 

Let us recall from Section IIV-AI that an (n, /.:) classical 
linear block code constructed over the code space C maps the 
information word x £ onto the corresponding codeword 
x £ F£. In the circuit-based representation, this encoding 
procedure can be encapsulated as follows: 

C = {x = (x : 0 n -k) V }, 


(64) 
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Fig. 13. Circuit representation of the inverse encoder eV 1 = (l : s). 


where V is an (n x n)-element invertible encoding matrix 
over F 2 and 0 n -k is an (n — /c)-bit vector initialized to 0. 
Furthermore, given the generator matrix G and the PCM H, 
the encoding matrix V may be specified as: 


V = 



and its inverse is given by: 


V- 1 = ( G - 1 H t ) . 


(65) 


( 66 ) 


The encoding matrix V specifies both the code space as well 
as the encoding operation, while its inverse V~ 1 specifies 
the error syndrome. More specifically, let y = x + e be the 
received codeword, where e is the n-bit error incurred during 
transmission. Then, passing the received codeword y through 
the inverse encoder V~ x yields: 

yV~ x = (x:s), (67) 


where x = x +1 for the logical error l £ F?j inflicted on the 
information word x and s £ F 2 _fe is the syndrome, which is 
equivalent to yH T . Eq. (l67l > may be further decomposed to: 

(x + e) V~ x = (x + l : s), 
xV~ l + eV~ l = {x : 0„_ fc ) + {l : s), (68) 


which is a linear superposition of the inverse of Eq. (l64t and 
eV~ x = (l : s ). Hence, the inverse encoder V~ 1 decomposes 
the channel error e into the logical error l and error syndrome 
s, which is also depicted in Fig. [13] 

Analogously to Eq. (El. the unitary encoding operation V 
of an [n, k] QSC, constructed over a code space C, which 
maps the information word (logical qubits) |V>) £ C 2 onto the 
codeword (physical qubits) | ip) £ C 2 ", may be mathematically 
encapsulated as follows: 


c = {|'0} = V(IV’) ® |0„_ fc »}, (69) 


where |0 n _fc) are (n — k) auxiliary qubits initialized to the 
state |0). The unitary encoder V of Eq. ( I69t i carries out an 
n-qubit Clifford transformation, which maps an n-qubit Pauli 
group Q n onto itself under conjugation II1051 , i.e. we have: 

vg n v^ = g n . ( 70 ) 


In other word, a Clifford operation preserves the elements of 
the Pauli group under conjugation such that for V £ Q n , 
VW' £ Q n . Furthermore, any Clifford unitary matrix is 
completely specified by a combination of Hadamard (H) gates. 


phase (S) gates and controlled-NOT (C-NOT) gates, which are 
defined as follows ED: 



C-NOT = 


1 

1 

/I 

0 

0 

VO 



0 0 0\ 
10 0 
0 0 1 
0 1 0 / 



(71) 


Hadamard gate preserves the elements of a single-qubit Pauli 
group 0 1 as follows: 

X -> HXH f = Z, 

Z HZH t = X, 

Y ->• HYH f = -Y, (72) 


while phase gate preserves them as: 

X ->• SXSt = Y, 

Z szs 1, = Z, 

Y -> SYS t = -X, (73) 


Since C-NOT is a 2-qubit gate, it acts on the elements of O 2 , 
transforming the standard basis of O 2 as given below: 


X ® I —>• X ® X, 

I ® X I ® X, 

Z ® I ->• Z ® I, 

I ® Z ->■ Z ( 8 ) Z. (74) 


Let us further emphasize on the significance of Clifford 
encoding operation. Since V belongs to the Clifford group, 
it presents the elements of the stabilizer group H under 
conjugation. If g[ is the ith stabilizer of the unencoded state 
\ip), then this may be proved as follows: 

\fj) ® |0 n _fc) = g\ (|V>) ® |0„_ fc )). (75) 

Encoding |i jj) with V yields: 

V (|^) ® |0 n _ fc )) = V {g[ (|tA) ® |0„_fc))), (76) 

which is equivalent to: 

V (|V>) ® |0„_ fc )) = V {g[V'V (|V») ® |0„- fe ))) , (77) 

since V'V = I n . Substituting Eq. (|69] > into Eq. (1771 ) gives: 

W) = ( Vg[V f ) m. (78) 

Hence, the encoded state \tp) is stabilized by g t = V//'VV From 
this it appears as if any arbitrary V (not necessarily Clifford) 
can be used to preserve the stabilizer subspace, which is not 
true. Since we assume that the stabilizer group jj is a subgroup 
of the Pauli group, we impose the additional constraint that V 
must yield the elements of Pauli group under conjugation as 
in Eq. <[70t . which is only true for Clifford operations. 

Furthermore, the Clifford encoding operation also presents 
the commutativity relation of stabilizers. Let g[ and gl be a 
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pair of unencoded stabilizers. Then the above statement can 
be proved as follows: 

9i9j = (Vs'Vt) = Vg'g'Vl (79) 

Since g[ and (/ J commute, we have: 

VgW^ = Vg' j g' i V t. (80) 

Using V^V = I„, gives: 

Vg'^Vg'Vl = g jgi . (81) 

Since the n-qubit Pauli group forms a basis for the (2 n x 2”)- 
element matrices of Eq. dZU, the Clifford encoder V, which 
acts on the 2 "-dimensional Hilbert space, can be completely 
defined by specifying its action under conjugation on the Pauli- 
X and Z operators acting on each of the n qubits, as seen in 
Eq. (l72l > to ( 1741 ). However, V and V', which differ only through 
a global phase such that V' = e je V, have the same impact 
under conjugation. Therefore, global phase has no physical 
significance in the context of Eq. d70l) and the n-qubit encoder 
V can be completely specified by its action on the binary 
equivalent of the Pauli operators. More specifically, for an n- 
qubit Clifford transformation, there is an equivalent 2 n x 2 n 
binary symplectic matrix V, which is given by: 


Zi,Xi ,..., Z . Xk are the encoded logical operators, which 
commute with the stabilizers g,. Logical operators merely map 
one codeword onto the other, without affecting the codespace 
C of the stabilizer code. It also has to be mentioned here that 
the stabilizer generators gi together with the encoded logical 
operations constitute the normalizer of the stabilizer code. The 
(2 n x 2 n)-element binary symplectic encoding matrix V is 
therefore given by: 


( z 1 \ 


/ z x \ 

Zk 


Zk 

Zk+i 


9i 

Z n 


Qn—k 

Xi 


Xi 

_Xfc 


x k 

Xfc+i 


tl 

\ / 


\k n —k) 


[VVV^] = [P]V = PV, (82) 

where [.] denotes the effective Pauli group G n such that P = 
[V] differs from V by a multiplicative constant, i.e. we have 
P = V /(± I ± *}, and the elements of G n are represented 
by 2n-tuple binary vectors based on the mapping given in 
Eq. @. As a consequence of this equivalence, any Clifford 
unitary can be efficiently simulated on a classical system as 
stated in the Gottesman-Knill theorem |106| . 

We next define V by specifying its action on the elements 
of the Pauli group Q n . More precisely, we consider 2 n n- 
qubit unencoded operators Zj .A,;...., Z n . X n , where Z, and 
Xi represents the Pauli Z and X operator, respectively, acting 
on the ith qubit and the identity I on all other qubits. The 
unecoded operators Zk+i ,..., Z n stabilizes the unencoded 
state of Eq. ( l69l >. i.e. ® |O ra _fc}), and are therefore 

called the unencoded stabilizer generators. On the other hand, 
Xk+i, ■ ■ ■, X n are the unencoded pure errors since X, anti¬ 
commutes with the corresponding unencoded stabilizer gen¬ 
erator Zi, yielding an error syndrome of 1. Furthermore, 
the unencoded logical operators acting on the information 
qubits are Z,. X,. .... Zk. Xf., which commute with the un¬ 
encoded stabilizers Zk+i,. ■ ■, Z n . The encoder V maps the 
unencoded operators Zi,Xi,..., Z n , X„ onto the encoded 
operators Zi, Xi ,..., Z n , X n , which may be represented as 
follows: 

X, = [VXiVi] = [X,] V, Zi = [VZiV t] = [Zi] V. (83) 

Since Clifford transformations do not perturb the commutativ¬ 
ity relation_of the operators, the resultant encoded stabilizers 
Zk+i, ■ ■ ■, Z n are equivalent to the stabilizers g, of Eq. ©, 
while Xfc+i,... , X n are the pure errors /:, of the resultant sta¬ 
bilizer code, which trigger a non-trivial syndrome. Moreover, 


where the Pauli Z and X operators are mapped onto the classi¬ 
cal bits using the Pauli-to-binary isomorphism of Section ITV-CI 
Analogously to the classical inverse encoder of Eq. (l67l >. the 
inverse encoder of a quantum code is the Hermitian conjugate 
)X. Let I?/)) = V\^) be the received codeword such that V is 
the n-qubit channel error. Then, passing the received codeword 
\ip) through the inverse encoder yields: 

V ] V\W) = v¥v(| ip) <8 |0 (n _*))) 

= {cm ® (5|0 (n _ fc) )), (85) 

where V^PV = (Cc/jS) and C (7 Gt- denotes the error imposed 
on the information word, while S G G n -k represents the error 
inflicted on the remaining (n — k) auxiliary qubits. In the 
equivalent binary representation, Eq. (l85l > may be modeled as 
follows: 

PV- 1 = {L:S), (86) 


where we have P = [P], L = [£] and S = [5], 

Let us now derive the encoding matrix V for the 3-qubit 
bit-flip repetition code, which has a binary PCM H given by: 


( 1 1 

0 

0 

0 

0 ^ 

V 1 0 

l 

0 

0 

0 ) 


(87) 


The corresponding encoding circuit is depicted in Fig. [14] Its 
unencoded operators are as follows: 
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(88) 
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Fig. 14. Encoding Circuit for 3-qubit bit-flip repetition code. 


A C-NOT gate is then applied to the second qubit, which 
is controlled by the first. As seen in Eq. (|74l). the C-NOT 
gate copies Pauli X operator forward from the control qubit 
to the target qubit, while Z is copied in the opposite direction. 
Therefore, we get: 
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(89) 


Another C-NOT gate is then applied to the third qubit, which 
is also controlled by the first, yielding: 
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= V. (90) 


As gleaned from Eq. ( l90t , the stabilizer generators of the 3- 
qubit bit-flip repetition code are gi = ZZI and g 2 = ZIZ. 
More explicitly, rows 2 and 3 of V constitutejhe PCM H of 
Eq . ®. The encoded logical operators are Z\ = ZII and 
X\ = XXX, which commute with the stabilizers <j\ and g 2 . 
Finally, the pure errors are t\ = IXI and t 2 = IIX, which 
anti-commute with gi and g 2 , respectively, yielding a non¬ 
trivial syndrome. 

Based on the above discussion, we now proceed to lay out 
the circuit-based model for a convolutional code, which is 
given in (331 . As discussed in Section IIV-DI convolutional 
codes are equivalent to linear block codes associated with 
semi-infinite block lengths. More specifically, as illustrated 
in Fig. [TT] the PCM H of an ( n,k,m ) convolutional code 
has a block-band structure, where the adjacent blocks have 
an overlap of m submatrices. Similarly, the encoder V of 
a classical convolutional code can be built from repeated 
applications of a linear invertible seed transformation U, which 
is an (n + m) x (n + m)-element encoding matrix, as shown 
in Fig. [15] The inverse encoder V~ x can be easily obtained 
by moving backwards in time, i.e. by reading Fig. [15] from 
right to left. Fet us further elaborate by stating that at time 
instant j, the seed transformation matrix U takes as its input 
the memory bits rrij- 1 £ F™, the logical bits lj £ ¥ 2 and the 
syndrome bits Sj £ F 2 _fc to generate the output bits ej £ F£ 



Fig. 15. Circuit representation of the encoder V of a convolutional code ED 


and the memory state rrij. More explicitly, we have: 

{rrij : ej) = (■ nrij _i : lj : Sj) U , (91) 

and the overall encoder is formulated as (33]: 

^ = U[l,...,n+m]U[n+l,...,2n+m] ■ • ■ U[(N—l)n+l,...,Nn+m]i 
N 

= 11 U[(j— l)n+l,...,jn-(-m]) (92) 

i=i 

where N denotes the length of the convolutional code 
and U [(j—i)n+i,...,m-(-m] acts on ( n + m ) bits, i.e. 
(rrij- 2 : lj -1 : Sj_i). For an [n, k,m\ quantum convolutional 
code, the seed transformation U is a 2 (n + m ) x 2 (n + m)- 
element symplectic matrix and Eq. ( |9TT > may be re-written as: 

(■ M j :Pj) = (Mj- 1 :Lj-.S j )U , (93) 

where M represents the memory state with an ?n-qubit Pauli 
operator. 

The aforementioned methodology conceived for construct¬ 
ing the circuit-based model of unassisted quantum codes may 
be readily extended to the class of entanglement-assisted 
codes (34]. The unitary encoding operation V of an [ n,k,c] 
EA-QSC, which acts only on the n transmitter qubits, may be 
mathematically modeled as follows: 

c = (|?) = V(\i/j) Tx <g> |0 q ) Tx (8) 1 0t) TxRx )}, (94) 

where the superscripts T\ and llx denote the transmitter’s 
and receiver’s qubits, respectively. Furthermore, |0 a ) Tx are a 
auxiliary qubits initialized to the state |0), where a = (n — 
k — c ), and \4>t) TxRx are the c entangled qubits. Analogously 
to Eq. (l85l) . the inverse encoder of an entanglement-assisted 
quantum code V gives: 

V f P|^) = V t PV(| ij}) Tx <8 |0 q ) Tx (8 | <Pt) TxRx ) 

= (C Tx \il>) Tx )® (5 Tx |0 a ) Tx ® (£ Tx \< t>t) TxRx ), 

(95) 

where C Tx £ Qk denotes the error imposed on the information 
word, while S Tx £ Q a represents the error inflicted on the 
transmitter’s a auxiliary qubits and £ lx £ Q c is the error 
corrupting the transmitter’s half of c ebits. The equivalent 
binary representation of Eq. ( |95| > is given by: 

PV- 1 =(L:S:E), 


(96) 
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where we have P = [' P Tx ], L = [C Tx ], S = [5 Tx ] and 
E = [E Tx ] . Similarly, Eq. J93b can be re-modeled as follows: 

(Mj : Pi) = {Mj-i : L, : Sj : E 3 ) U. (97) 


B. System Model: Concatenated Quantum Codes 

Fig. [T6l shows the general schematic of a quantum com¬ 
munication system relying on a pair of concatenated quantum 
stabilizer codes. In this contribution, both the inner as well 
as the outer codes are assumed to be convolutional codes. 
Furthermore, analogously to the classical concatenated codes, 
the inner code must be recursive, while both the inner as well 
the outer code must be non-catastrophic. Having a recursive 
nature of the inner code is essential for the sake of ensur¬ 
ing that the resultant families of codes have an unbounded 
minimum distance. On the other hand, the non-catastrophic 
nature of both the inner and the outer codes guarantees that 
a decoding convergence to an infinitesimally low error rate 
is achieved. It was found in |39l , |98l that QCCs cannot 
be simultaneously recursive and non-catastrophic. In order to 
overcome this problem, Wilde et al. l34l . Il36l proposed to 
employ entanglement-assisted inner codes, which are recursive 
as well as non-catastrophic. Therefore, the inner code should 
be an entanglement-assisted recursive and non-catastrophic 
code, while the outer code can be either an unassisted or an 
entanglement-assisted non-catastrophic code. 

At the transmitter, the intended quantum information \ipi) 
is encoded by an [rii,fci] outer encoder Vi using (m — k\) 
auxiliary qubits, which are initialized to the state |0), as 
depicted in Eq. ( l69l >. The encoded qubits |i/q) are passed 
through a quantum interleaver (7r). The resultant permuted 
qubits |i/) 2 ) are fed to an [ 712 ,^ 2 ] inner encoder V 2 , which 
encodes them into the codewords |r£o} using ( 712 —£ 2 ) auxiliary 
qubits initialized to the state |QrT The n-qubit codewords 
|i/> 2 }> where we have n = 711712 , are then serially transmitted 
over a quantum depolarizing channel, which imposes an 71 - 
tuple error V 2 £ Gn on the transmitted codewords. 

At the receiver, the received codeword |'02) = V 2 IV’ 2 ) is 
passed through the inverse encoder v|, which yields the cor¬ 
rupted information word of the inner encoder £ 2 ^ 2 ) and the 
associated (712 — fc 2 )-qubit syndrome <S 2 | 0 („ 2 _fc 2 )) as depicted 
previously in Eq. (f85l> . where C 2 denotes the error imposed 
on the logical qubits of the inner encoder, while S> represents 
the error inflicted on the remaining (712 — fe) qubits. The 
corrupted logical qubits of the inner encoder are de-interleaved, 
resulting in Vi\4>i), which is then passed through the inverse 
outer encoder V-j. This gives the corrupted information word 
of the outer encoder C\\ipi) and the associated (711 — £:i)-qubit 
syndrome <Si|0 (ni _ fel) ). 

The next step is to estimate the error C\ for the sake of 
ensuring that the original logical qubit | ipi) can be restored by 
applying the recovery operation TZ. For estimating C 1 , both the 

--Please note that this is a general schematic. The inner code can be either 
an un-assisted or an entanglement-assisted code. However, it is advisable to 
use an entanglement-assisted inner code for the sake of ensuring an unbounded 
minimum distance of the resultant concatenated code. 


syndromes 6 > 2 | 0 („ 2 _fc 2 )} and <Si | 0 ( ni _ fcl )) are fed to the inner 
and outer Soft-In Soft-Out (SISO) decoders ll27ll . respectively, 
which engage in iterative decoding ll33l . f34ll in order to yield 
the estimated error C\. The corresponding block is marked as 
‘MAP Decoder’ in Fig. [16] Here, P"(.), P?(.) and P°(.) denote 
the a-priori , extrinsic and a-posteriori probabilities IZ7ft related 
to the zth decoder. Based on this notation, the turbo decoding 
process can be summarized as follows: 

• The inner SISO decoder of Fig. [16] uses the chan¬ 
nel information P c h(7 ? 2), the a-priori information 
gleaned from the outer decoder P'jfZA) (initialized 
to be equiprobable for the first iteration) and the 
syndrome S 2 to compute the extrinsic information 
P 2 (£ 2 ). For a coded sequence of length N, we 
have V 2 = [7\i,^ 2 , 2 , ■ • •, V 2 ,u ■ • •, P2,n\> where 

V2,t = [Pi v ^2 t> ■ • • i p2t\- The channel information 
Pch (/ 2 .t) is computed assuming that each qubit is 
independently transmitted over a quantum depolarizing 
channel having a depolarizing probability of p, whose 
channel transition probabilities are given by lf33l : 


{pit) 


1 -P, if Pl t = I 

p/3, ifp* it e{x,z,Y}. 


(98) 


• P 2 (P 2 ) is passed through the quantum de-interleaver 
(7r _1 ) of Fig. [16] to generate the a-priori information 
for the outer decoder Pi('Pi). 

• Based on both the a-priori information Pi('Pi) and on 
the syndrome <Si, the outer SISO decoder of Fig. [16] 
computes both the a-posteriori information P) J (jC i ) and 
the extrinsic information Pj {V \). 

• Pi (Pi) is then interleaved to obtain P 2 (£ 2 ), which is fed 
back to the inner SISO decoder of Fig. [16] This iterative 
procedure continues, until either convergence is achieved 
or the maximum affordable number of iterations is 
reached. 

• Finally, a qubit-based MAP decision is made for deter¬ 
mining the most likely error coset C\. It must be men¬ 
tioned here that both the inner and outer SISO decoders 
employ the degenerate decoding approach of ll33l . which 
aims for finding the ‘most likely error coset’ rather than 
the ‘most likely error’ acting on the logical qubits £,, 
as we will discuss in the next section. 


C. Degenerate Iterative Decoding 

As discussed in Section IIV-BI quantum codes exhibit the 
intrinsic property of degeneracy, which is also obvious from 
Eq. (l85l >. More explicitly, we have: 

5|0„_ fc )=5 1 |0)(8)---®5 n _fc|0). (99) 

Since, we have Si £ {I, X,Y, Z}, we can re-write Eq. ( l99b 
as follows |33| : 


5|0„_fe) = e|si) 0 • • • ® | Sn — k) 1 


( 100 ) 
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Fig. 16. System Model: Quantum communication system relying on concatenated quantum stabilizer codes. P“(.), P®(.) and P?(.) denote the a-priori, extrinsic 
and a-posteriori probabilities related to the ith decoder. 


where e £ {±l,±i}, and: 

Si = 0 if Si = I or Si = Z, 

Si = 1 otherwise. (101) 

For example, if .Si = Y and 5, = I for i / 1, since Y = ?'XZ, 
we get S\0 n - k ) = *|1) ® |0) <g> • • ■ <g) |0). 

Observing the (n — k) syndrome qubits of Eq. < 1 100b col¬ 
lapses them to the classical syndrome s = {si, ..., s n -fe}, 
which is equivalent to the symplectic product of P and H, i.e. 
s = (P * //,), More precisely, the syndrome sequence 

|0„_fc) is invariant to the Z-component of S since Z|0) = |0). 
Let S be the effective 2(n— fc)-bit error on the syndrome, which 
may be decomposed as S = S x + S z , where S x and S z are the 
X and Z components of .S', respectively. Then s only reveals 
S x . Hence, two distinct error sequences P = (L : S x + S Z )V 
and P' = (L : S x + S' Z )V, which only differ in the Z- 
component of S, yield the same syndrome s. Furthermore, it 
must be noted that both P and P' have the same logical error 
L. Therefore, P and P' differ only by the stabilizer group 
and are known as degenerate errors, which do not have to be 
distinguished, since they can be corrected by the same recovery 
operation L~ l . 

Recall that a classical syndrome-based MAP decoder aims 
for finding the most likely error for a given syndrome, which 
may be modeled as: 

L{S) = argmax i P(L|S l ), (102) 

where P(L|5) denotes the probability of experiencing the 
logical error L imposed on the transmitted qubits, given that 
the syndrome of the received qubits is S. By contrast, quantum 
codes employ degenerate decoding, which aims for finding the 
most likely error coset C(L, S x ) associated with the observed 
syndrome S x . The coset C(L, S x ) is defined as [[33): 

C(L, S x ) = {P=(L:S X + S' 2 )!/} VS Z £ {I, Z} n ~ k . 

(103) 



P e (L) 

P e (P) 


Fig. 17. General schematic of a SISO decoder. P“(.), P e (.) and P°(.) denote 
the a-priori, extrinsic and a-posteriori probabilities. 


Therefore, a degenerate MAP decoder yields: 

L{S X ) = argmax i P(L|5 ,a: ), (104) 

where we have: 

P(L\S X ) = P(L\(S X +S Z )). (105) 

S z G{I,Z} n ~ k 

The MAP decoder of Fig. [16] consists of two serially con¬ 
catenated SISO decoders, which employ the aforementioned 
degenerate decoding approach. Fig. [17] shows the general 
schematic of a SISO decoder, where the Pauli operators V, 
C and S are replaced by the effective operators P, L and 
S x , respectively. The SISO decoder of Fig. [T7] yields the a- 
posteriori information pertaining to L and P based on the 
classic forward-backward recursive coefficients a and /3, as 
follows [|33l : 

• For a coded sequence of duration N, let 
P = [P 1 ,P 2 ,...,P t ,...,P N ] and L = 

[L 1 ,L 2 ,...,L t ,...,L N ], where P t £ G n and 
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L t G Gk • More explicitly, P t = [P*, P t 2 ,..., P t n ] and 

• Let us decompose the seed transformation as U = ( Um '■ 
Up ), where Um is the binary matrix formed by the first 
2m columns of U, while Up is the binary matrix formed 
by the last 2 n columns of U. Therefore, we have: 

M t = (Mt-i : L t : S t ) U M , (106) 

P t = (M t - 1 :L t :St)U P . (107) 

• Let at (Mt) be the forward recursive coefficient, which 
is defined as follows: 

at (Mt) — P (M t \S< t ) ) 

cx Y, P ° = A ) P ° ( P *) “t-i M - 008 ) 

where £| t = (S?) Q<j<t , p G G m , A G G k and 
a G G n -k, while a ~=~ a x + a z , having a x = Sf. 
Furthermore, we have P t = (/i : A : <j)Up and 
M t = (p : A : a) U M . 

• Let [3t (M t ) be the backward recursive coefficient, which 
is defined as: 

Pt (M t ) = P(M t \S> t ), 

cxJ2 pa ( p * = A ) pa ( p t+ 1) A+1 (Aft+1), 

A,cr 

(109) 

where -Pt+i = (M t : A : a)U P and 

M t+ i = (M t : A : ct)Um- 

• Finally, we have the a-posteriori probabilities P°(L t ) 
and P°(P t ), which are given by: 

P°(L t )±P(L t \S x ), 

oc J> a (L t )P a (P t H-i (»& (M t ), (110) 

/A,cr 

P°(Ps) — P(P:|<S' X ), 

a P Q (P t )P a (L t = A)a t _r ( M )/3 t (Af t ), 

/i,,A ,cr 

(ill) 

where S’ 21 = (S'f) 0 < t < JV , P t = (p : L t : a)U P and 
M t = (p: L t : a) U M . 

• The marginalized probabilities P°(Lj), for j G {0. /:: — 
1}, and P °(Pt), for j G {0 ,n — 1}, are then computed 
from P°(Lj) and P°(P/), respectively. The a-priori in¬ 
formation is then removed in order to yield the extrinsic 
probabilities |[34l . i.e we have: 

ln[P e (Lj)] = ln[P °(Ll)} - ln[P“(^)], (112) 

ln[P e (P/)] = ln[P°(P/)] - ln[P a (P/)]- (113) 

It has to be mentioned here that the property of degeneracy 
is only an attribute of auxiliary qubits and the ebits of an 
entanglement-assisted code do not contribute to it. This is 


because both X as well as Z errors acting on the transmitter’s 
half of ebits give distinct results when measured in the Bell 
basis, i.e. £ Tx \<fit) TxRx gi ves four distinct Bell states for 
£j x G {I, X, Z,Y}. Consequently, the degeneracy is a 
function of a and reduces to zero for a = 0. 

VI. EXIT-Chart Aided Code Design 

EXIT charts 1271 . |32l , |107| are capable of visualizing 
the convergence behaviour of iterative decoding schemes by 
exploiting the input/output relations of the constituent decoders 
in terms of their average Mutual Information (MI) character¬ 
istics. The EXIT chart analysis not only allows us to dispense 
with the time-consuming Monte-Carlo simulations, but also 
facilitates the design of capacity approaching codes without 
resorting to the tedious analysis of their distance spectra. 
Therefore, they have been extensively employed for designing 
near-capacity classical codes HSU-HID- Let us recall that 
the EXIT chart of a serially concatenated scheme visualizes 
the exchange of four MI terms, i.e. average a-priori MI of the 
outer decoder 1\, average a-priori MI of the inner decoder 
l\, average extrinsic MI of the outer decoder I [ p: , and average 
extrinsic MI of the inner decoder l\. More specifically, l\ 
and constitute the EXIT curve of the outer decoder, while 
l\ and l\ yield the EXIT curve of the inner decoder. The 
MI transfer characteristics of both the decoders are plotted in 
the same graph, with the x and y axes of the outer decoder 
swapped. The resultant EXIT chart quantifies the improvement 
in the mutual information as the iterations proceed, which can 
be viewed as a stair-case-shaped decoding trajectory. An open 
tunnel between the two EXIT curves ensures that the decoding 
trajectory reaches the (l,y) point of perfect convergence. 

In our prior work (35], we extended the application of 
EXIT charts to the quantum domain by appropriately adapting 
the conventional non-binary EXIT chart generation technique 
for the quantum syndrome decoding approach. Recall from 
Section IIV-CI that a quantum code is equivalent to a classical 
code. More specifically, the decoding of a quantum code is 
essentially carried out with the aid of the equivalent classical 
code by exploiting the additional property of degeneracy, as 
discussed in Section IV-CI Quantum codes employ syndrome 
decoding, which yields information about the error-sequence 
rather than about the information-sequence or coded qubits, 
hence avoiding the observation of the latter sequences, which 
would collapse them back to the classical domain. Since a 
quantum code has an equivalent classical representation and 
the depolarizing channel is analogous to a Binary Symmetric 
Channel (BSC), we employ the EXIT chart technique to design 
hashing bound approaching concatenated quantum codes. The 
major difference between the EXIT charts conceived for the 
classical and quantum domains is that while the former models 
the a-priori information concerning the input bits of the inner 
encoder (and similarly the output bits of the outer encoder), 
the latter models the a-priori information concerning the 
corresponding error-sequence, i.e. the error-sequence related 
to the input qubits of the inner encoder L 2 (and similarly the 
error-sequence related to the output qubits of the outer encoder 

Pi). 
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Fig. 18. System model for generating the EXIT chart of the inner 
decoder ED 



Fig. 19. System model for generating the EXIT chart of the outer 
decoder ED 


Similar to the classical EXIT charts, it is assumed that the 
interleaver length is sufficiently high to ensure that ll27l . If32l : 

• the a-priori values are fairly uncorrelated; and 

• the a-priori information has a Gaussian distribution. 

Fig. 1 1 81 shows the system model used for generating the EXIT 
chart of the inner decoder. Here, a quantum depolarizing 
channel having a depolarizing probability of p generates the 
error sequence P 2 , which is passed through the inverse inner 
encoder . This yields both the error imposed on the logical 
qubits L‘z and the syndrome Sf. The a-priori channel block 
then models the a-priori information P“(Z,2) such that the 
average MI between the actual error L 2 and the a-priori 
probabilities is given by I A {L 2 ) (27) , (32) , |107| . More 

explicitly, we have I A {L 2 ) = /[T 2 , where I denotes 

the average MI function. Moreover, the ith and (N + i)th 
bits of the effective error vector can be visualized as 
4-ary symbols. Consequently, similar to classical non-binary 
EXIT charts (97), 01 121 , the a-priori information is modeled 
using an independent Gaussian distribution with a mean of 
zero and variance of a\, assuming that the X and Z errors 
constituting the 4-ary symbols are independent Based on 
the channel information P ch ( P>) , on the syndrome Sf and on 
the a-priori information, the inner SISO decoder generates the 
extrinsic information P 2(^2) by using the degenerate decoding 
approach of Section IV-CI Finally, the extrinsic average MI 
-Te(P 2 ) = /[T2, P2(^2)] between L 2 and is computed. 

Since the equivalent classical capacity of a quantum channel 
is given by the capacity achievable over each half of the 4-ary 
symmetric channel, is the normalized MI of the 4-ary 

symbols, which can be computed based on (97), 111 131 as: 


EP2(^ (m) )l0g 2 P^(4 (m) )]V 

,m=0 J / 

( 1 } 4 ) 

where E is the expectation (or time average) operator and L^' 
is the m lh hypothetical error imposed on the logical qubits. 
More explicitly, since the error on each qubit is represented 
by an equivalent pair of classical bits, L'i}' 1 is a 4-ary classical 


Ie(L 2 ) = - 2 


23 Under the idealized asymptotic conditions of having an infinite-length 
interleaver, I a (£ 2 ) may be accurately modeled by the Gaussian distribution. 
As and when shorter interleavers are used, the Gaussian assumption becomes 
less accurate, hence in practice a histogram-based approximation may be relied 
upon. 


symbol associated with m G {0, 3}. The process is repeated for 
a range of /a(T 2 ) G [0,1] values for the sake of obtaining the 
extrinsic information transfer characteristics at the depolarizing 
probability p. The resultant inner EXIT transfer function T 2 of 
the specific inner decoder may be defined as follows: 

I E (L 2 )=T 2 [I A (L 2 ),p\, (115) 

which is a function of the channel’s depolarizing probability 
P- 

The system model used for generating the EXIT chart of 
the outer decoder is depicted in Fig. [19] As inferred from 
Fig. [19] the EXIT curve of the outer decoder is independent 
of the channel’s output information. The a-priori information 
is generated by the a-priori channel based on Pi (error on 
the physical qubits of the second decoder) and I A (Pi), which 
is the average MI between Pi and P“(Pi). Furthermore, as 
for the inner decoder. Pi is passed through the inverse outer 
encoder Vff 1 to compute Sf, which is fed to the outer SISO 
decoder to yield the extrinsic information P®(Pi). The average 
MI between Pi and P®(Pi) is then calculated using Eq. (II 14k 
The resultant EXIT chart is characterized by the following MI 
transfer function: 


Ie(Pi) = Ti[I a (Pi)], (116) 

where Ti is the outer EXIT transfer function, which is depen¬ 
dent on the specific outer decoder, but it is independent of the 
depolarizing probability p. 

Finally, the MI transfer characteristics of both decoders 
characterized by Eq. (II 15b and Eq. ( 1116b are plotted in the 
same graph, with the x and y axes of the outer decoder 
swapped. For the sake of approaching the achievable capacity 
of Fig. [3] our EXIT-chart aided design aims for creating a 
narrow, but marginally open tunnel between the EXIT curves 
of the inner and outer decoders at the highest possible depo¬ 
larizing probability (analogous to the lowest possible SNR for 
a classical channel). For a given noise limit p* and the desired 
code parameters, this may be achieved in two steps. We first 
find that specific inner code, which yields the largest area under 
its EXIT-curve at the noise limit p*. Once the optimal inner 
code is selected, we find the optimal outer code, whose EXIT- 
curve gives the best match with the chosen inner code. The 
narrower the tunnel-area between the inner and outer decoder’s 
EXIT curve, the lower is the deviation from the achievable 
capacity, which may be quantified using Eq. ©. 
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VII. A Key to Hashing Bound: Quantum 
Irregular Convolutional Codes 

In this section, we exploit the EXIT-chart aided design 
criterion of Section [VI] to design concatenated codes, which 
operate arbitrarily close to the hashing bound. Here, we assume 
that we already have the optimal inner code. More explicitly, 
our design objective is to find the optimal outer code C 
having a coding rate Rq, which gives the best match with the 
given inner code, i.e. whose EXIT curve yields a marginally 
open tunnel with the given inner decoder’s EXIT curve at 
a depolarizing probability close to the hashing bound. For 
the sake of achieving this objective, a feasible design option 
could be to create the outer EXIT curves of all the possible 
convolutional codes to find the optimal code C, which gives the 
best match, as we did in our prior work f35l . To circumvent 
this exhaustive code search, in this contribution we propose to 
invoke Quantum Irregular Convolutional Codes (QIRCCs) for 
achieving EXIT-curve matching. 

Similar to the classical Irregular Convolutional Code (IRCC) 
of QH, our proposed QIRCC employs a family of Q sub¬ 
codes C q , q £ {1,2,..., Q}, for constructing the target code 
C. Due to its inherent flexibility, the resultant QIRCC provides 
a better EXIT-curve match than any single code, when used as 
the outer component in the concatenated structure of Fig. IT6l 
The q th subcode has a coding rate of r q and it encodes a 
specifically designed fraction of the original information qubits 
to g q N encoded qubits. Here, N is the total length of the 
coded frame. More specifically, for a Q-subcode IRCC, o q is 
the q th IRCC weighting coefficient satisfying the following 
constraints 033, EH: 

Q Q 

^Qq = 1 , R Q = '^2,eq' r q . Qq S [0,1], Vq , (117) 

9=1 9=1 

which can be conveniently represented in the following matrix 
form: 

(n k ::: / e ) <« ^> T = (i) 

r e = R . (118) 

Hence, as shown in Fig. [20] the input stream is partitioned 
into Q sub-frameQ which are assembled back into a single 
iV-qubit stream after encoding. 

In the context of classical IRCCs, the subcodes C q are 
constructed from a mother code 111 141 . Ill 151 . More specifically, 
high-rate subcodes are obtained by puncturing the mother 
code, while the lower rates are obtained by adding more 
generators. However, unlike classical codes, puncturing is 
not easy to implement for quantum codes, since the resul¬ 
tant punctured code must satisfy the symplectic criterion, as 
in HI 161 . In this context, in order to design the constituent 
subcodes of our proposed QIRCC, we selected 5 strong 
randomly-constructed memory-3 quantum convolutional codes 

24 This is only true if all subcodes are active. If o q = 0 for the q th subcode, 
then C q is not activated. Therefore, the input stream is only divided among 
the active subcodes. 



Fig. 20. Structure of a Q-subcode QIRCC encoder. 

with quantum code rates {1/4,1/3,1/2, 2/3,3/4}, which met 
the non-catastrophic criterion of (33]. More explicitly, for 
the sake of achieving a random construction for the Clifford 
encoder specifying the quantum convolutional code, we used 
the classical random walk algorithm over the (n + m)-qubit 
Clifford group as in HI 17H . The seed transformations of the 
resultant subcodes having rates {1/4,1/3,1/2, 2/3, 3/4} are 
given below: 

U\ = {9600,691,11713,4863,1013,6907,1125,828,10372, 
6337,5590,11024,12339,3439}, 

U 2 = {3968,1463,2596,3451,1134,3474,657,686,3113, 
1866,2608,2570}, 

U 3 = {848,1000,930,278,611,263, 744,260,356,880}, 

U 4 = {529,807,253,1950,3979,2794,956,1892,3359,2127, 
3812,1580}, 

U 5 = {62,6173,4409,12688, 7654,10804,1763,15590,6304, 
3120,2349,1470,9063,4020}. (119) 

The EXIT curves of these QIRCC subcodes are shown in 
Fig. [2X1 whereby the memory-3 subcodes of Eq. ( 1119b are 
indicated by solid lines. Furthermore, in order to facilitate 
accurate EXIT curve matching with a sufficiently versatile 
and diverse set of inner EXIT functions, we also selected 
5 weak randomly-constructed memory-1 subcodes for the 
same range of coding rates, i.e. {1/4,1/3,1/2, 2/3, 3/4}. The 
corresponding seed transformations are as follows: 

U 6 = {475,194,526,422,417,988,426,611,831,84}, 

U 7 = {26,147,149,99,112,184,64,139}, 

U 8 = {37,55,58,35,57,54}, 

U 9 = {57,248,99,226,37,93,244,54}, 

U w = {469,634,146,70,186,969,387,398,807,452}, 

( 120 ) 

and their EXIT curves are plotted in Fig. [2T] with the aid 
of dotted lines. It must be mentioned here that the range of 
coding rates chosen for the QIRCC subcodes can be expanded 
such that the EXIT curves cover a larger portion of the EXIT 
plot, which further improves curve matching. However, this 
increases the encoding and decoding complexity. 

Based on our proposed QIRCC, relying on the 10 subcodes 
specified by Eq. ( II 1 9b and (1 1 20b . the input bit stream is 
divided into 10 fractions corresponding to the 10 different- 
rate subcodes. The specific optimum fractions to be encoded 
by these codes are found by dynamic programming. More 
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qircc-subcodeslO.gle 



Fig. 21. Outer EXIT curves (inverted) of our QIRCC subcodes having code 
rates {1/4,1/3,1/2, 2/3, 3/4} for both memory-3 as well as memory-1. 


specifically, since the QCCs belong to the class of linear 
codes, the EXIT curves of the 10 subcodes, given in Fig. [2TJ 
are superimposed onto each other after weighting by the 
appropriate fraction-based weighting coefficients, which are 
determined by minimizing the area of the open EXIT-tunnel. 
To elaborate a little further, the transfer function of the QIRCC 
is given by the weighted sum of each subcode’s transfer 
function as shown below: 


Q 

Ie(Pi) = Ti[/ a (Pi)] = ^ g q T* [I A (Pi)], (121) 

9=1 


where we have: 

/T 2 [i 1 ,p}\ 

T2[l2,P] 


\T 2 [*n,p]/ 



Tl l [h] 
Tf 1 N 


T^iiA 

T?~\h] 


(124) 


W _1 [*n] Tf 1 ^] ... T?-\i n ]j 


Here, N denotes the number of sample points such that i € 
{ii,i 2 , ■ ■ ■ ■ in} and it is assumed that N > Q. Furthermore, 
the error should be greater than zero for the sake of ensuring 
an open tunnel, i.e. we have: 


e(i) > 0, Vi g [0,1], (125) 


The resultant cost function, i.e. sum of the square of the errors, 
is given by EE): 

J{Qi,-->,Qq)=[ e(i) 2 di, (126) 

Jo 

which may also be written as: 

J(Q) = e T e. (127) 


The overall process may be encapsulated as follows: 

Qo P t = ar g min J{q), (128) 

e 

subject to Eq. (1117l > and (1 1 25b . which is a convex optimization 
problem. The unconstrained optimal solution for Eq. ( 1 1 28b 
is found iteratively using steepest descent approach with a 
gradient of dj(g)/dg = 2e, which is then projected onto 
the constraints defined in Eq. (II1 7b and (11251 ). Further details 
of this optimization algorithm can be found in 111 141 . HI 151 . 


VIII. Results and Discussions 


where T/ [Ia (Pl)] is the transfer function of the q th subcode. 
For a given inner EXIT curve and outer code rate Rq, 
we employ the curve matching algorithm of EH. EH 
for optimizing the weighting coefficients g of our proposed 
QIRCC such that the square of the error between the inner and 
inverted outer EXIT curves is minimized subject to Eq. dim . 
More explicitly, the error function may be modeled as: 

e(i) =T 2 [i,p\-T^ 1 [i], (122) 

where p = (p* — e) given that p* is the noise limit defined 
by the hashing bound and e is an arbitrarily small number. 
The corresponding matrix-based notation may be formulated 

as EH- EE): 

(123) 


For the sake of demonstrating the curve matching capability 
of our proposed QIRCC, we designed a rate-1/9 concatenated 
code relying on the rate-1/3 entanglement-assisted inner code 
of f34j , ll36l . namely “PTOIREA”, with our proposed QIRCC 
as the outer code. Since the entanglement consumption rate of 
“PTOIREA” is 2/3, the resultant code has an entanglement 
consumption rate of 6/9, for which the corresponding noise 
limit is p* = 0.3779 according to Eq. © E). Furthermore, 
since we intend to design a rate-1/9 system with a rate-1/3 
inner code, we have Rq = 1/3. Hence, for a target coding 
rate of 1/3, we used the optimization algorithm discussed in 
Section IVIII for the sake of finding the optimum weighting 
coefficients of Eq. (1128b at the highest possible depolarizing 
probability p = p*—e. It was found that we only need to invoke 
two subcodes out of the 10 possible subcodes, based on g = [0 
0 0 0 0.168 0.832 0 0 0 0] T , for attaining a marginally open 
tunnel, which occurs at p = 0.345, as shown in Fig. [22] Hence, 
the resultant code has a convergence threshold of p = 0.345, 


e = b — A^, 
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Depolarizing Probability 

Fig. 23. WER performance curves with increasing iteration number for an 
interleaver length of 3, 000 qubits. Rate-1/9 concatenated code, relying on 
PTOIREA as the inner code and the proposed QIRCC as the outer code, was 
used. 


Fig. 22. EXIT curves of the concatenated rate-1/9 system, with PTOIREA 
as the inner code and QIRCC as the outer, at p = 0.345 and p = 0.34. 


which is only [10 x log 10 ( o 0 3 3 7 4 7 5 9 )] = 0.4 dB from the noise 
limit of 0.3779. Fig. [22] also shows two decoding trajectories 
at p = 0.34 for a 30,000 qubit long interleaver. As gleaned 
from the figure, the decoding trajectories closely follow the 
EXIT curves reaching the (1,1) point of perfect convergence. 

The corresponding Word Error Rate (WER) performance 
curves recorded for our QIRCC-based optimized design using 
a 3,000 qubit long interleaver are seen in Fig. [23] where 
the WER is reduced upon increasing the number of itera¬ 
tions. More explicitly, our code converges to a low WER 
for p < 0.345. Thus, this convergence threshold matches the 
one predicted using EXIT charts in Fig. [22] More explicitly, 
since the EXIT chart tunnel closes for p > 0.345, the system 
fails to converge, if the depolarizing probability is increased 
beyond 0.345. Hence, the performance does not improve upon 
increasing the number of iterations if the depolarizing proba¬ 
bility exceeds the threshold. By contrast, when the depolarizing 
probability is below the threshold, the WER improves at 
each successive iteration. It should also be noted that the 
performance improves with diminishing returns at a higher 
number of iterations. 

Fig. l24l compares our QIRCC-based optimized design with 
the rate-1/9 “PTOIREA-PTOIR” configuration of l34| , which 
is labeled “A” in the figure. An interleaver length of 3000 
qubits was used. For the “PTOIREA-PTOIR” configuration, 
the turbo cliff region emerges around 0.31, which is within 0.9 
dB of the noise limit. Therefore, our QIRCC-based design out¬ 
performs the “PTOIREA-PTOIR” configuration of l34l . More 
specifically, the “PTOIREA-PTOIR” configuration yields a 
WER of 10” 3 at p = 0.29, while our design gives a 



Depolarizing Probability 

Fig. 24. Comparison of WER performance of our QIRCC-based optimized 
rate-1/9 QTC with the PTOIREA-PTOIR configuration of f34l (labeled “A”) 
and the exhaustive-search based optimized QTC of l35l (labeled “B”) for an 
interleaver length of 3, 000 qubits and a maximum of 15 iterations. 


WER of 10 3 at p = 0.322. Hence, our optimized design 
outperforms the ‘PTOIREA-PTOIR” configuration by about 
[10 x log 10 (£jj|)] = 0.5 dB at a WER of 1CT 3 . It must 
be mentioned here that the “PTOIREA-PTOIR” configuration 
may have a lower error floor than our design, yet our design 
exhibits a better performance in the turbo cliff region. We 
further compare our QIRCC-based optimized design with the 
exhaustive-search based optimized turbo code of |[35l , which 
is labeled “B” in Fig. [24] Both code designs have similar 
convergence threshold. However, our QIRCC-based design has 
a much lower error rate, resulting in a lower error floor as 
gleaned from Fig. [24] 
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IX. Conclusions and Design Guidelines 

Powerful QECCs are required for stabilizing and protecting 
the fragile constituent qubits of quantum computation as well 
as communication systems against the undesirable decoher¬ 
ence. In line with the developments in the field of classical 
channel coding theory, this may be achieved by exploiting 
concatenated codes designs, which invoke iterative decoding. 
Therefore, in this paper we have laid out a slow-paced tu¬ 
torial for designing hashing bound approaching concatenated 
quantum codes using EXIT charts. To bridge the gap between 
the quantum and classical channel coding theory, we have 
provided insights into the transition from the classical to the 
quantum code design. More specifically, with the help of toy 
examples, we have illustrated that quantum block codes as 
well as convolutional codes may be constructed from arbitrary 
classical linear codes. We then move onto the construction 
of concatenated quantum codes, focusing specifically on the 
circuit-based structure of the constituent encoders and their 
equivalent classical representation as well as the degenerate 
iterative decoding. Finally, we have detailed the procedure for 
generating EXIT charts for quantum codes and the principles 
of EXIT-chart aided design. Our design guidelines may be 
summarized as follows: 

• As discussed in the context of our design objectives in 
Section m we commence our design by determining 
the noise limit p* for the desired code parameters, i.e 
the coding rate and the entanglement consumption rate 
of the resultant concatenated quantum code, which was 
introduced in Section 11 

• We then proceed with the selection of the inner stabilizer 
code of Fig. [T6] which has to be both recursive as well 
as non-catastrophic, as argued in Section IV-BI Since 
the unassisted quantum codes cannot be simultaneously 
both recursive as well as non-catastrophic, we employ an 
entanglement-assisted code. Furthermore, the EA inner 
code of Fig. [16] may be either derived from the family of 
known classical codes, as discussed in Section [IV] or it 
may be constructed using random Clifford operations, 
which were discussed in Section IV-AI At this point, 
the EXIT curves of Section [VI] may be invoked for the 
sake of finding that specific inner code, which yields the 
largest area under its EXIT-curve at the noise limit p*. 

• Finally, we find the optimal non-catastrophic outer code 
of Fig. [16] which gives the best EXIT-curve match 
with that of the chosen inner code. In this context, 
our EXIT-chart aided design of Section [VI] aims for 
creating a narrow, but marginally open tunnel between 
the EXIT curves of the inner and outer decoders at the 
highest possible depolarizing probability. The narrower 
the tunnel-area, the lower is the deviation from the 
hashing bound, which may be quantified using Eq. ©. 

Recall that the desired code structure may also be optimized 
on the basis of a range of conflicting design challenges, which 
were illustrated in Fig. [4] 

Furthermore, for the sake of facilitating the hashing bound 
approaching code design, we have proposed the structure of 


QIRCC, which constitutes the outer component of a con¬ 
catenated quantum code. The proposed QIRCC allows us to 
dispense with the exhaustive code-search methods, since it 
can be dynamically adapted to match any given inner code 
using EXIT charts. We have constructed a 10-subcode QIRCC 
and used it as an outer code in concatenation with a non- 
catastrophic and recursive inner convolutional code of (34), 
(361 , In contrast to the concatenated codes of l34l . whose per¬ 
formance is within 0.9 dB of the hashing bound, our QIRCC- 
based optimized design operates within 0.4 dB of the noise 
limit. Furthermore, at a WER of 10 :i , our design outperforms 
the design of l34l by around 0.5 dB. Our proposed design also 
yields lower error rate as compared to the exhaustive-search 
based optimized design of lf35ll . 
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